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The unit of nonlinear interaction in Navier-Stokes and the Magnetohydrodynamic (MHD) equa- 
tions is a wavenumber triad (k,p,q) satisfying k + p + q = 0. The expression for the combined 
energy transfer from two of these wavenumbers to the third wavenumber is known. In this paper we 
introduce the idea of an effective energy transfer between a pair of modes through the mediation of 
the third mode and then find an expression for it. In fluid turbulence, energy transfer takes place 
between a pair of velocity modes, whereas in MHD turbulence energy transfer takes places between 
(1) a pair of velocity modes, (2) a pair of magnetic modes, and (3) between a velocity and a magnetic 
mode in a triad. In this paper we have obtained the expression for each of these transfers. We also 
show how the effective mode-to-mode energy transfer rate can be utilised to study energy cascades 
and shell-to-shell energy transfer rates. 
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I. INTRODUCTION 



In fluid and MHD turbulence, eddies of various sizes interact amongst themselves; energy is exchanged among 
them in this process. These interactions arise due to the nonlinearity present in these systems. The fundamental 
interactions involve wavenumber triads (k,p,q) with k + p + q = 0. 

The combined energy transfer computation to a mode from the other two modes of a triad has generally been 
considered to be fundamental. This formalism has played an important role in furthering our understanding of locality 
in fluid turbulence Jl]-[l| , and also in analysing subgrid scale eddy viscosity • Energy interactions are considered 
to be local if the triads with k ~ p ~ q contribute most to the energy transfer to the wavenumber k. Numerical 
simulations of Domaradzki and Rogallo showed that the contribution of the nonlocal triads [k ~ p > > q) to the 
interactions is not small, but in such triads the dominant exchange of energy occurs between wavenumbers of similar 
magnitudes, i.e., between k and p. Domaradzki et al. || computed the sub-grid scale viscosity by summing up the 
energy lost by a wavenumber k < k c (cut off wavenumber) by triad interactions with atleast one wavenumber larger 
than k c . Batchelor M first suggested a formula for the energy transfer in fluid turbulence between two regions in 
Fourier space, and Domaradzki and Rogallo jlj calculated shell-to-shell energy transfer rates. Taking energy transfer 



between shells as an example, we will show in Section [II, that these formalisms have certain limitations. 

In this paper we present a scheme to calculate the energy transfer rates between two modes in a triad interaction. 
Using our scheme we can calculate the energy transfer rates between two shells, cascade rates, and also study locality of 
energy transfer. The formulae for the combined energy transfer rate were used to study energy transfer rates in MHD 
turbulence [p|-p"2"|. Ishizawa and Hattori |[2) numerically computed the kinetic energy flux to the velocity modes 
outside the sphere and the magnetic energy fluxes to the magnetic modes outside the sphere. They also computed 
the kinetic energy fluxes due to transfer of energy to all magnetic modes and a similar magnetic energy flux due to 
transfer of energy to all velocity modes 12 1. Pouquet et al. ||, Pouquet JhJ, and Ishizawa and Hattori used 
the combined energy transfer formalism in an Eddy damped quasi-normal Markovian (EDQNM) closure calculation 
to study non-local energy transfer. However, several energy transfers cannot be calculated using these formalisms. 
Our formalism of "mode-to-mode" transfer enables us to calculate all these energy transfer rates between the spheres 
and the shells. 

This paper is organized into the following sections. In Section || we will state the known results regarding the 
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combined energy transfer between velocity modes of a triad. In Section III we will present a discussion on the formulae 



which are used for studying energy transfer between shells in fluid turbulence. In Section IV we will formulate a way 
of studying energy transfer between a pair of modes within a triad, with the third mode of the triad meditating 
the transfer. This new method is used to re-define shell-to-shell transfer in Section The known results of energy 



transfer in a MHD triad are discussed in Section VI and the "mode-to-mode" transfers in MHD in Section VII. The 



results of Section VII are used in Section VIII to define shell-to-shell transfer rates and the energy cascade rates in 



MHD turbulence. The conclusions of this paper follows in Section IX 



II. ENERGY TRANSFER IN A TRIAD: KNOWN RESULTS 



The Navier Stokes equation in real space are written as 



du 



^ + (u.V)u = -Vp + iA7 2 u, 

at 



(1) 



where u is the velocity field, and v is the fluid kinematic viscosity. In Fourier space, the kinetic energy equations for 
a Fourier mode is 



dE u (k) 
dt 



2vk 2 E"(k)= ^""(k|p,q), 



(2) 



k+p+q=0 



where E u (k) = |u(k)| /2 is the kinetic energy in Fourier space for mode k. The nonlinear term ^""(klp, q) is 



^"(klp.q) = -Re (i[k.u(q)][u(k).u(p)] +i[k.u(p)][u(k).u(q)]) 



(3) 



This term represents the combined transfer of kinetic energy from mode p and q to mode k |13| . Note that the 
wavenumber triad k, p, and q should satisfy the condition k + p + q = 0. This nonlinear term satisfies the following 
conservation property 



^""(klp, q) + S uu (p\k, q) + 5 M "(q|k, p) = 0. 



(4) 



While the quantity S""(k|p, q) represents the nonlinear energy transfer from the two modes p and q to mode k, it 



would be useful to know the exact energy transfer between any two modes, say from p to k. In Section IV we will 
explore this idea. 
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III. SHELL-TO-SHELL ENERGY TRANSFER IN FLUID TURBULENCE 



Using the combined energy transfer S uu (k|p, q), Domaradzki and Rogallo |jj have discussed energy transfer between 
two shells. They interpret the quantity 



as the rate of energy transfer from shell n to shell m [^jjj]- Note that k-sum is over shell m, p-sum over shell n, 
and q = k — p. However, Domaradzki and Rogallo themselves points out that it may not be entirely correct to 
interpret the formula (^) as the shell-to-shell energy transfer. The reason for this is as follows: 

In the energy transfer between two shells m and n, two types of wavenumber triads are involved (see Fig. 1 shown 
below). In a triad of type I, the wavenumbers p, q are located in one shell, and k is located in the other. In a type 
II triad, wavenumber k is in one shell, p or q in the other shell, and the third wavenumber is located outside the 
two shells. In Eq. (^) the summation is carried over both kinds of triads. However, it is easy to see that the energy 
transfer from shell n to shell m takes place through both the k-p and k-q legs of triad I but only through the k-p 
leg of triad II. Hence Batchelor's and Domaradzki's formalism do not yield correct shell-to-shell energy transfers (as 
was pointed out by Batchelor and Domaradzki themselves). 



The nonlinear interaction in the Navier-Stokes and the MHD equations are intrinsically three- mode interactions. 
The expression for the energy transfer to one mode of the triad from the other two was presented in Section ||. 
In this section we shall explore the possibility of obtaining an expression for the energy transfer between any two 
modes within a triad — we will call it the mode-to-mode transfer. We emphasize that this approach is still within the 
framework of the triad interaction — that is, the triad is still the fundamental interaction of which the mode-to-modc 
transfer is a part. The energy transfer between two modes of a triad by the mediation of the third mode is sought 
here. 

We shall first consider only the Navier-Stokes equation. In a later section, the discussion will be generalised to the 
MHD equation. 




(5) 



kem pen 



IV. MODE-TO-MODE ENERGY TRANSFER IN A TRIAD 
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A. Definition of Mode-to-Mode transfer in a triad 



Consider a triad k, p, q. Let the quantity ^. u "(k|p|q) denote the energy transferred from mode p to mode k with 
mode q playing the role of a mediator (see Fig. 2). We wish to obtain an expression for lfi au . 
The ljt uu 's should satisfy the following relationships : 

1. The sum of ^. uu (k|p|q) and ^ u "(k|q|p), which represent energy transfer from mode p to mode k and from 
mode q to mode k, respectively, should be equal to the total energy transferred to mode k from modes p and 
q, i.e., S uu (k\p,q) [see Eq. (§)]. Thus, 

# uu (k|p|q) + ^"(klqlp) = S«"(k|p, q), (6) 



# u "(p|k|q) + ^""(Plqlk) = S uu (p|k, q), (7) 

^"(qlklp) + ^""(qlPlk) = 5""(q|k, p). (8) 

2. The energy transferred from mode p to mode k, i.e., ^" u (k|p|q), will be equal and opposite to the energy 
transferred from mode k to mode p i.e., ^. uu (p|k|q). Thus, 

r"(k|p|q)+r"(p|k|q)=0, (9) 



r u (k|q|p) + r u (q|k| P ) = o, (10) 
r u (p|q|k) + r m (q|p|k) = o. (ii) 

These are six equations with six unknowns. However, the value of the determinant formed from the Eqs. @-(|ll]) 
is zero. Therefore we cannot find unique T/t uu, s given just these equations. However, it is reasonable to expect that 
there is a definite amount of energy transfer from one mode to the other mode, say from mode k to mode p given 
u(k), u(p), u(q). Since Eqs. do not yield a unique 1fc au , we need to use constraints based on invariance, 

symmetries, etc. to get a definite T/t uw using Eqs. (||)-(|l|). 



5 



B. Solutions of equations of mode-to-mode transfer 

One solution of Eqs. (|)-@ is 

r"(k|p|q) ee -Re (i[k.u(q)] [u(k).u(p)]) , (12) 
From the definition of ^ u "(k|p|q), it directly follows that $ uu, s satisfy the following conditions. 

r"(k|p|q) + r "(k|q|p) = 5"«(k|p, q), (13) 

^"(Plklq) +r"(p|q|k) = S™(p|k,q), (14) 

r u (q|k| P ) + ^""(qlPlk) = S«"(q|k, p). (15) 

Using the triad relationship k + p + q = 0, and the incompressibility constraint [k.u(k) =0], it can be seen that 
$ uu, s, etc. also satisfy the following conditions 

r"(k| P |q)+r tl (p|k|q) = 0, (16) 

r u (k|q| P )+r"(q|k|p)-0, (17) 

r u (p|q|k)+r"(q|p|k) = 0. (18) 



Comparing Eqs. (jig) -(jig) with Eqs. (||)-(11), it is clear that the set of $ uu, s is one instance of the T/t uu, s : i.e., 
^. Mtl (k|p|q) = ^""(k|p|q). However, this is not a unique solution. 

If another solution ^ UIl (k|p|q) differs from ^ u "(k|p|q) by an arbitrary function X A , i.e., ^t""(k|p|q) = ^" u (k|p|q)+ 
Xa, then by inspection we can easily see that every possible solution of Eqs. (](j)-(jll]) must be of the form 

r"(k|p|q)=r tl (k|p|q)+X A (19) 
r"(k|q|p)=r"(k|q|p)-X A (20) 

r"(p|k|q)=r"(p|k|q)-X A (21) 
^ tl (p|q|k)=^ u (p| q |k)+X A (22) 



r u (q|k|p) = r tl (q|k| P ) + x A 



(23) 



r u (q|p|k) = r u (q|p|k)-x A 



(24) 



Note that Xa can depend upon the wavenumber triad k, p, q, and the Fourier components u(k), u(p), u(q). A a 
needs to be determined from other symmetry and invariance arguments to uniquely fix T/t wu . 

In Appendix A we have attempted to determine Xa- We construct Xa by observing that it can depend on k, p, 
q, u(k), u(p), u(q); and it must satisfy rotational invariance, galilean invariance, and it should be finite. In our 
procedure we write down a general scalar function of k, p, q, u(k), u(p), and u(q) with undetermined coefficients. 
The coefficients can be written as a series in the scalars of the type k.p, k.u(p), u(k).u(p) with unknown constants. 
By imposing invariance, and finiteness constraints we have shown that, upto the linear order in these scalars, the 
coefficients must vanish, giving A a = 0. Unfortunately, we have not been able to determine the higher order 
coefficients. Hence, from the analysis of Appendix A we cannot show that Aa will definitely vanish. However, if 
Aa = we get a simple relationship, 



giving the energy transfer between a pair of modes in a triad with the third mode mediating the transfer, as simply 
$ uu . But if A a ^ then the energy transfer picture is more complex and the energy transfer between a pair of modes 
will be given by Eqs. (|l9])-(^4|) where Aa is unknown. However, even if Xa is unknown we can obtain very relevant 
information from the $ uu, s as is shown below. 

We will now give a physical interpretation to the two parts of Tjt uu, s, i.e., $ MU 's and Aa. We now see from Fig. 3 
that 5? u "(k|p|q) + Aa gets transferred from p to k, y? u "(q|k|p) + Xa gets transferred from k to q, j#" u (p|q|k) + Xa 
gets transferred from q to p. The quantity Aa flows along p — > k — > q — > p, circulating around the entire triad 
without changing the energy of any of the modes. Therefore we will call it the Circulating transfer. Of the total 
energy transfer between two modes, $ uu + Xa, only $ uu can bring about a change in modal energy. Aa transferred 
from mode p to mode k is transferred back to mode p via mode q, i.e., the mode p transfers Xa directly to mode k, 
and mode p transfers Aa back to k indirectly through mode q. Thus the energy that is effectively transferred from 
mode p to mode k is just y?""(k|p|q). Therefore 5# u "(k|p|q) can be be termed as the effective mode-to-mode transfer 
from mode p to mode k. 



Tji 



UU 



,11 II 



(25) 
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In summary, we have attempted to obtain an expression for the energy transfer rate Tjt between two modes in 
a triad. To leading order in a series expansion, we can show that 1jt uu — $ wu . But the ambiguity in TjL u remains 
because of the lack of a complete proof. However, the important conclusion is that the mode-to-mode energy transfer 
can be expressed as T/l uu = $ uu + Xa where $ uu is the 'effective mode-to-mode transfer' and is the 'circulating 
transfer'. In the next section we will further develop the notion of effective transfers by applying it to energy transfer 
between shells. 



V. SHELL-TO-SHELL ENERGY TRANSFER AND CASCADE RATES USING "MODE-TO-MODE" 

FORMALISM 

In Section III we had pointed out the problem in treating the expression in Eq. (^|) as the shell-to-shell transfer. 
This expression, commonly taken to be the shell-to-shell transfer, actually gives the energy transferred to shell m from 
both shell n and the modes q of type II triads (see Fig. 1). In this section we will make use of the idea of "circulating" 
and "effective mode-to-mode transfer" to redefine shell-to-shell transfer. 

Consider energy transfer between shells m and n in Fig. 4 . The mode k is in shell m, the mode p is in shell n, and 
q could be inside or outside the shells. In terms of the mode-to-mode transfer "J?. uu (k|p|q) from mode p to mode k, 
the energy transfer from shell n to shell m can be defined as 

T r n -EE^ u ( k ip|q) ( 26 ) 

kern pen 

where the k-sum is over shell m, p-sum is over shell n, and k+p+q^O. The quantity T/t wu can be written as a sum 
of an effective transfer .# utl (k|p|q) and a circulating transfer X<\- We know from the last section that the circulating 
transfer does not contribute to the energy change of modes. From Fig. 4 we can see that X& flows from shell m to 
shell n and then flows back to m indirectly through the mode q. Therefore the effective energy transfer from shell m 
to shell n is just ^ uu (k|p|q) summed over all k in shell m and all p in shell n, i.e., 

kem pen 

where ^ u "(k|p|q) is the effective mode-to-mode transfer. 

If Xa = then the effective energy transfer between two shells given by this equation will be the same as the same 
as the total energy transfer between them given by Eq. (pq). However, we have argued earlier that the circulating 
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transfer does not result in the energy change of any of the modes. Hence, from the physical point of view the effective 
shell-to-shell transfer is a very relevant quantity to study. In a subsequent paper we have studied the effective 
shell-to-shell transfers in numerical simulations. 

A. Energy cascade rates 

The kinetic energy cascade rate (II) (or flux) in fluid turbulence is defined as the rate of loss of kinetic energy by 
a sphere in k-space to the modes outside. In literature, the energy flux in fluid turbulence has been computed using 
5""(k|p, q) |r^ , [l5|] . Kraichnan M] shows by direct integration of Eq. (||) that the kinetic energy lost from a sphere 
of radius K by nonlinear convective transfers is 

n W = - E E ^(k|p,q). (28) 

|k|<K|p|>A' 

Although the energy cascade rate in fluid turbulence can be found by the above formula, the mode-to-mode approach 
provides a more natural way of looking at the energy transfers. In later sections we will obtain expressions for 
mode-to-mode transfer in MHD turbulence and define various fluxes which are not accessible to the conventional 
approach. Here, we will obtain an expression for the flux in terms of the effective mode-to-mode energy transfer. 
Since ^t u "(k|p|q) represents energy transfer from p to k with q as the mediator, we may alternatively write the 
energy loss from a sphere as 

n W = - E E r"(k| P |q). (29) 
|k|<K |p|>K 

where mode k is inside the sphere, p is outside the sphere, and q = k — p. The mode-to-mode transfer ^t u "(k|p|q) 
consists of a circulating part and an effective part. From Fig. 5 we see that the net amount of circulating transfer 
leaving the sphere is zero. Thus the circulating transfer makes no contribution to the the energy flux from the sphere. 
Therefore, it can be removed from Eq. ( p9| ) and the resultant expression for the flux can be written as 

n (*0 = - E E r u ( k ip|q)- (so) 

|k|<RT|p|>if 

with a summation over the effective mode-to-mode transfer. The expressions in Eqs. ( p8| ) and (|3^) are equivalent as 
is formally proved in Appendix B. 
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VI. ENERGY TRANSFER IN A MHD TRIAD: KNOWN RESULTS 



The MHD equations in real space are written as 

du 



+ (u.V)u = -Vp+(b.V)b + ^V 2 u, (31) 
at 



and 



<9b 

— + (u.V)b= (b.V)u + M V 2 b, (32) 

where u and b are the velocity and magnetic fields respectively, and v and /i are the fluid kinematic viscosity and 
magnetic diffusivity, respectively. In Fourier space, the kinetic energy and magnetic energy evolution equations for a 
Fourier mode are 

r)E u (]<\ 1 1 

—±l + 2vk 2 E u (k) = ]T -S""(k|p,q)+ J2 2^ (k|P ' q) ' (33) 

k+p+q=0 k+p+q=0 

^M+2M 2 ^(k)= £ ±S»(k|p,q) + ^(k|p,q), (34) 

' k+p+q=0 k+p+q=0 

where E u (k) = |u(k)| /2 is the kinetic energy, and E b (k) = |b(k)| 2 /2 is the magnetic energy. The four nonlinear 
terms 5""(k|p,q), 5" b (k|p,q) S^^k^.q) and 5 hu (k|p,q) are 

S™(k|p,q) = -i?e(i[k.u(q)][u(k).u(p)] + i[k.u(p)][u(k).u(q)]) , (35) 

5 bb (k|p,q) = -i?e(i[k.u(q)][b(k).b(p)] + z[k.u(p)] [b(k).b(q)]) , (36) 

5" b (k|p,q) =i?e(i[k.b(q)][u(k).b(p)] + i[k.b(p)][u(k).b(q)]) , (37) 

5 b "(k|p,q) =i2e(i[k.b(q)][b(k).u(p)] +i[k.b(p)][b(k).u(q)]). (38) 



These terms are conventionally taken to represent the nonlinear transfer from modes p and q to mode k |l^Jl_3j of a 
triad such that k + p + q = 0. The term S uu (k\p, q) represents the net transfer of kinetic energy from modes p and 
q to mode k. Likewise the term 5'" & (k|p, q) is the net magnetic energy transferred from modes p and q to the kinetic 
energy in mode k, whereas S bu (k\p, q) is the net kinetic energy transferred from modes p and q to the magnetic 
energy in mode k. The term S' bf, (k|p, q) represents the transfer of magnetic energy from modes p and q to mode k. 
All these transfer terms are represented in the Fig. 6. 
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Stanisic |16j showed that the nonlinear terms satisfy the following detailed conservation properties: 



5""(k|p, q) + S uu (p|k, q) + 5 uu (q|k, p) = 0, 



(39) 



S bb (k|p,q) + S bb (p|k,q) + S bb (q|k,p 



■) = o, 



(40) 



and 



S" b (k|p, q) + S" b (p|k, q) + S" b (q|k, p) + S bu {k\p, q) + S bu {p\k, q) + S bu (q|k, p) = 0. 



(41) 



The Eq. (|39j) implies that kinetic energy is transferred conservatively between the velocity modes of a wavenumber 
triad, and the Eq. ( [40| ) implies that magnetic energy is also transferred conservatively between the magnetic modes 
of a wavenumber triad. The Eq. (p~[| ) implies that the cross transfers of kinetic and magnetic energy, S ub (k\p, q) and 
S bu (k\p, q), within a triad are also energy conserving. 

The quantities 5 u "(k|p,q), 5 ui> (k|p,q), 5 6b (k|p,q), and 5 bu (k|p,q) represent the nonlinear energy transfer from 
the two modes p and q to mode k. It would be useful to know the energy transfer between any two modes, say from 
p to k. This issue will be explored in the next section. 



In Section ^ we had presented a new approach to describe the energy transfer between two velocity modes of a 
triad mediated by the third velocity mode, in the Navier-Stokes equation. In the same spirit we will now attempt 
to find the mode-to-mode energy transfers in MHD. In MHD turbulence there will be three kinds of mode-to-mode 
transfer within a wavenumber triad k, p, q: kinetic energy transfer from u(p) to u(k); magnetic energy transfer from 
b(p) to b(k); and transfer of kinetic energy from u(p) to magnetic energy in b(k). We denote these three transfers 
by ^ Mtl (k|p|q), ^?. bb (k|p|q), and ^. bu (k|p|q) respectively, where the index q indicates that the energy transfer between 
modes p and k is mediated by the mode q. Since the nonlinear interactions fundamentally involve three-modes, the 
energy transfer between a pair of modes should depend on the third mode. The mode-to-mode transfers 7fc uu , ljl bb , 
and 1fi bu are schematically illustrated in Fig. 7. 

In this section we will obtain a description for each of these transfers. 



VII. "MODE-TO-MODE" ENERGY TRANSFERS IN MHD EQUATIONS 
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A. Velocity mode to velocity mode energy transfers 

In Section [i^ we discussed the mode-to-mode transfer, 1JL UX1 , between velocity modes in fluid flows. In this section 
we will find T/t uu for MHD flows. The transfer of kinetic energy between the velocity modes is brought about by the 
term (u.V)u, both in the Navier Stokes and MHD equations. Hence, the expression for the combined kinetic energy 
transfer to a mode from the other two modes of the triad is also same for the two. i.e., the combined transfer to u(k) 
from u(p) and u(q) is given by 5 uu (k|p,q) [see Eqs. (|) and ©]. Consequently, TjL uu 's for MHD will satisfy the 
constraints given in Eqs. (|^)-(|ll|) for the corresponding l^ uu 's for fluids. As a result, ^t""(k|p,q) in MHD can be 
expressed as a sum of a circulating transfer and the effective transfer ^ utl (k|p|q) given by Eq. (|lj), i.e., 

r"(k|p|q)=r tl (k|p|q)+X A (42) 



The arguments are the same as those presented in Section IV B for fluid turbulence. 



B. Magnetic mode to Magnetic mode energy transfers 

Now we consider the magnetic energy transfer from mode b(p) to b(k) in the triad (k,p,q) [see Fig. 7]. This 
transfer which is denoted by ^. bb (k|p|q), should satisfy the following relationships : 

1. The sum of the mode-to-mode energy transfers from b(p) to b(k) and from b(q) to b(k) should be equal to 



S bb (k\p, q) [Eq. ( pq ) in Section VI], which is the combined energy transfer to b(k) from b(p) and b(q). Thus, 

^ fcb (k|p|q) + ^ bb (k|q|p) = S bb (k\p, q), (43) 

^ bb (p|k|q) + ^ bb (p|q|k) = S bb (p\k, q), (44) 

^ bb (q|k|p) + # bb (q|p|k) = S" b (q|k, p). (45) 



2. the energy transfer from b(k) to b(p), ^?. bb (k|p|q), and the transfer from b(p) to b(k), ^. bb (p|k|q), should be 



equal but opposite in sign. That is, 



^ bb (k|p|q) + ^ bb (p|k|q)=0, (46) 



^ bb (k|q|p)+^ bb (q|k|p)=0, (47) 
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^ bb ( P |q|k)+^ bb (q| P |k)=0. 



(48) 



Once again, the above equations cannot uniquely determine ^?. bb (k| P |q) since the value of the determinant formed 
from these equations is zero. While discussing the mode-to-mode energy transfers in Navier Stokes equation in 
Section IV we got the same result for ^. utl (k|p|q). Following the reasoning in Section |y|, we can extract physically 
meaningful information from Eqs. ([43[)-(f48|). 



The combined energy transfer to b(k) from b(p) and b(q) is given by Eq. (36) of Section VI. We denote the first 



term on the right hand side of that equation by ^ bb (k|p|q), i.e 



Ak|p|q) = -ite(i[k.u(q)][b(k).b(p)]), 



(49) 



By replacing $ bb (k\p\q) for ^t 0D (k|p|q) in Eqs. (|43|)-(|48|) we find that ^ D0 (k|p|q) is a solution of these equations, i.e., 



ibb. 



bb i 



^ b (k| P |q)+^°(k|q| P )=5 bb (k| P ,q), 



?bb 



(50) 



^ bb (p|k|q)+^ ( P |q|k)=5 6f, (p|k,q), 



,66 



bb, 



(51) 



^ bb (q|k| P )+^ D (q| P |k) = ^(q|k, P 



ybb 



(52) 



and 



^ bb (k|p|q)+^ bb ( P |k|q) = 0, 



(53) 



^ bb (k|q| P )+^ t '(q|k| P ) = 0, 



bb, 



(54) 



^ bb ( P |q|k)+/ b (q| P |k) = 0. 



(55) 



Thus, ># bb (k| P |q) is a solution of the equations ( [43|) — (|4£|). By inspection it can be seen that all solutions of the 
equations can be expressed as 



^ bb (k|p|q)= Ak|p|q)+y A , 



(56) 



:Ak|q| P )= Ak|q|p)-y A , 



^ bb ( P |k|q)=^°( P |k|q)-y A , 



66 / 



(57) 
(58) 
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^ fcb (p|q|k)= Ap|q|k) + F A , 



(59) 



^ fcb (q|k|pH Aq|k|p) + y A , 



(60) 



^ b (q|p|k)=^ b (q|p|k)-y A , 



(61) 



where Ya is an arbitrary scalar function dependent on the wavenumber triad k, p, q and the Fourier components 
u(k), u(p), u(q), b(k), b(p), b(q) at those wavenumbers. To determine the unique form of this function we 
need to impose additional constraints, similar to those imposed on X/\ in Appendix A for determining mode-to-mode 



energy transfer between velocity modes in fluid turbulence (see Section IVB). The arguments for determination of 
Ya is given in Appendix B. 

The solutions in Eqs. (p6|)-(|6l|) have been schematically illustrated in Fig. 8. We see from the figure that Ya is 
transferred from b(p) — > b(k) — > b(q) — > and back to b(p), i.e., it circulates around the triad without causing any 
change in modal energy — it is a thus a circulating transfer. The magnetic energy effectively transferred from b(p) 
to b(k) is just <# bb (k|p|q)[],i.e., 



^ b / (k| P |q)=^ bb (k|p|q). 



(62) 



The effective mode-to-mode transfer ^ bb (k|p|q) from b(p) to b(k) is mediated by the velocity mode u(q). 

In the next section we will use T^Jfj to calculate effective energy transfer rate between the magnetic modes inside 
two shells and the effective cascade rate of magnetic energy. 



C. Velocity mode to Magnetic mode energy transfers 



In Section VII A and Section VII B we discussed mode-to-mode energy transfer between a pair of velocity modes 
and between a pair of magnetic modes in a wavenumber triad, respectively. We now consider the energy transfer 
between u(p) and b(k), ^. Ilb (k|p|q), within the triad (k, p, q) illustrated in Fig. 7 . We will follow the same 



sequence of steps as in the two sections VII A and VII B 



T^, ub 's will satisfy the following relationships: 



Hhe idea of circulating transfer and effective mode-to-mode transfer in a triad was introduced in Section 
discussion can be found in that section. 
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IVB 



— a detailed 



1. Since ^" 6 (k|p|q) and ^" 6 (k|q|p) are the mode-to-mode energy transfers from b(p) to u(k) and from b(q) to 
u(k), the sum of the two should be equal to S"" b (k|p, q), the combined energy transfer to u(k) from b(p) and 
b(q). Therefore, we get 

r b (k|p|q) + r 6 (k|q|p) = S ub (k\p, q), (63) 
^(Plklq) + ^(Pklk) = S ub (p\k, q), (64) 
r b (q|k| P ) + ^(qlPlk) - S" b (q|k, p). (65) 



^ 6 "(k|p|q) + ^"(k|q|p) = S fe "(k| P , q), (66) 



r"(p|k|q) + ^"(p|q|k) = S fc "(p|k, q), (67) 



^(q|k|p) +^(q|p|k) = S btt (q|k,p). (68) 



2. "^" b (k|p|q) denotes the energy transfer from b(p) to u(k). ^. hu (p|k|q) is physically the same transfer but seen 
as a transfer from u(k) to b(p). Therefore, 

r b (k|p|q)+^ btl ( P |k|q)=0, (69) 



r b (k|q|p) + ^ btl (q|k|p)=0, (70) 



r b (p|q|k) + ^ fc "(q|p|k)-o. (71) 



^ fc "(k|p|q) + r b (p|k|q)=0, (72) 



^(k|q|p)+r b (q|k| P )-0, (73) 



^"(p|q|k)+r b (q|p|k)-0. (74) 
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The solutions of these equations are not unique and the expression for the energy transfer between modes cannot 
be obtained from the above equations alone. As in sections |y| and VII A, we will now explore the solutions of the 
above equations. 

We define the following quantities : 

r & (k|p|q) = -Re (i[k.b(q)] [u(k).b(p)]) , (75) 



^(k|p|q) = -Re ( i [k.b(q)][b(k).u(p)]) , (76) 
Replacing $ ub and $ bu for lfc Llh and 7fc bu respectively in Eqs. (^3|)-(0) we find that the are a solution to those 



equations, i.e. 



and 



r 6 (k|p|q) + r b (k|q|p) = S llb (k|p, q), (77) 
^(Plklq) + ^(Plqlk) = S ub (p\k, q), (78) 



r b (q|k| P ) + ^(qlPlk) = 5" b (q|k, p). (79) 



£ 6u (k|p|q) + ^ b "(k|q|p) = S bu (k\p, q), (80) 



^(Plklq) + ^ b "(p|q|k) = S bu (p\k } q), (81) 



^"(q|k|p) + ^"(q|p|k) = S bu (q\k, p), (82) 



r b (k| P |q)+^ btl (p|k|q)-0, (83) 



r b (k|q| P )+^ bu (q|k|p) = 0, (84) 



r b (p|q|k)+^ btl (q| P |k) = 0, (85) 



^"(k|p|q)+r b (p|k|q) = 0, (86) 
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^"(k|q|p)+r b (q|k|p)-0, (87) 



^(p|q|k) + r b (q|p|k) = 0. (88) 

The $ ub 's are just a single instance of the the T/t ub, s. It can be seen by inspection that all solutions can be expressed 
in the form : 

^ U (k|p|q)=^"(k|p|q) + Z A , (89) 
^"(k|q|p)=^"(k|q|p)-Z A , (90) 



# 6 >|k| q )=£ 6 >|k|q)-Z A , (91) 



^ b "(p|q|k)^ b "(p|q|k) + Z A , (92) 



^"(q|k| P )=^"( q |k|p) + Z A , (93) 



^ b "(q|p|k) ^ b "(q|p|k) - Z A , (94) 



^ b (k|p|q)=^(k|p|q) + Z A , (95) 



# u6 (k|q|p)=£ tt6 (k|q|p)-Z A , (96) 



^ b (p|k|q)=^(p|k|q)-^A, (97) 



r"(p|q|k) = r 6 (p|q|k) + ^A, (98) 



# u6 (q|k|p) = r b (q|k| P ) + Z A , (99) 



r h (q|p|k) = r h (q|p|k)-^A, (ioo) 

where Z\ is an arbitrary function dependent on k, p, q, u(k), u(p), u(q), b(k), b(p), b(q). These solutions are 
pictorially represented in Fig. 9 below. 
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We are already familiar with such solutions from the earlier discussions. From Fig. 9 we see that Z\ transfers 
energy from u(p) — > b(k) — > u(q) — > b(p) — > u(k) — > b(q) — > and back to u(p) without resulting in a change in modal 
energy. Hence, following the discussions in Section VII A and Section VII B| , we can interpret the quantity as a 



circulating transfer. The <fi bu can be interpreted as the effective mode-to-mode transfers. For example, <# b "(k|p|q) is 
the effective transfer from u(p) to b(k),i.e, 

^ / (k|p|q)=^ u (k|p| q ), (101) 

which is mediated by the magnetic mode b(q). 

In the next section, we will use the effective mode-to-mode transfer to define, (a) shell-to-shell transfers between 
the velocity and the magnetic modes, and (b) cascade rates between the velocity and the magnetic modes. 



VIII. SHELL-TO-SHELL ENERGY TRANSFER AND CASCADE RATES IN MHD TURBULENCE 

In this section, we will extend the formulation of Section |y| to define shell-to-shell transfer rates and cascade rates 
in MHD turbulence. 

We shall use the term w-shell and M-sphere respectively to denote a shell and a sphere in fc-space containing velocity 
modes, and 6-shell and 6-sphere for the corresponding magnetic modes. 



A. Shell-to-Shell energy transfer rates 



In Section [v| we had defined the effective shell-to-shell transfer rates between two it-shells in Fourier space. The 
expression for the effective mode-to-mode transfer between two velocity modes is the same for both Fluid and MHD 



flows [see Section IV B and VII A | . Hence, the effective shell-to-shell transfer rates between two it-shells in MHD will 
also be defined by Eq. ( p7| ) of Section ^, i.e., 



kem pen 



(102) 



Similarly, the quantity ^ b6 (k|p|q) gives the energy transfer rate from b(p) to b(k), mediated by q. Hence the 



energy transfer rate from n th 6-shell to the m th 6-shell can be obtained by summing ^t bb (k|p|q) over every p in the 



the n th 6-shell and over every k in the vn th 6-shell. From Section VII B we know that 7jt bb — ^ bb (k|p|q) + where 



5# bb (k|p|q) is the effective mode-to-mode transfer from b(p) to b(k), and Ya is a circulating transfer. Thus, is 
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transferred from b(p) in shell n to b(k) in shell m and then back to b(p) via the mode b(q). Therefore, the effective 
shell-to-shell transfer transfer rate from the n th to the m th 6-shell is 



i£ fe n = EEAk| P |q). ( 103 ) 



kem pen 



We can now also dehnc shell-to-shell transfers rates between a u-shcll and a 6-shell. Tjt u (k|p|q) is the energy transfer 
rate from u(p) to b(k) mediated by mode q. Hence the energy transfer rate from the the n th u-sheil to the m th 6-shell 
can be calculated by summing ^. b "(k|p|q) over every p in the n th it-shell and over every k in m th 6-shell. We have 
shown in Section VII C that ^?. b "(k|p|q) = fi bu (k\p\q) + where ^ b "(k|p|q) is the effective mode-to-mode transfer 
from the mode u(p) to the mode b(k) in the triad k, p, q, and is the circulating transfer which is transferred along 
u(p) — > b(k) — > u(q) — > b(p) — > u(k) — > b(q) — > u(p). Hence, Za is transferred from the mode u(p) in shell n to 
b(k) in shell m and then flows back to u(p) via the modes u(q) — » b(p) — » u(k) — ► b(q) — > u(p) — » u(p). Therefore, 
we interpret the quantity obtained by summing 5# bu (k|p|q) over every p in the n th u-shell and every k in the m th 
6-shell as the effective shell-to-shell transfer. That is, 

^n-EE^WPl^) ( 104 ) 
kern pen 

is the effective transfer from the the n th w-shell to the m th 6-shell. The energy transfer rate from the n th 6-shell to 
the m th u-shell, T^ n = -T% n . 

In Appendix A-C, we have argued that the circulating transfer A"a = 0, Ya = 0, and — to the first order in 
the expansion coefficients. If — 0, Ya = 0, Z& — 0, then the effective shell-to-shell transfers will be the same as 
the total shell-to-shell transfers. However, we had pointed out is Section [v| that the circulating transfer, if nonzero, 
does not contribute to a change in energy of any shell, and hence from the physical point of view the effective shell- 
to-shell transfers is of greater significance than the total shell-to-shell transfers. We have numerically computed the 



effective shell-to-shell transfers given by Eqs. (102)- (104) in numerical simulations of 2-D MHD turbulence and gained 
important insights into the energy transfer process. The results of the simulations will be discussed in a companion 
paper @. 



B. Energy cascade rates 



There are various types of cascade rates (energy fluxes) in MHD turbulence. We have schematically shown these 
transfers in Fig. 10. In this section we will derive the formulae to calculate these cascade rates within the framework 
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of effective transfers. 



In Section V A, we derived an expression for kinetic energy flux in terms of the mode-to-mode transfer [Eq. (|26|)]. 
We showed that the circulating transfer Xa does not contribute to the flux. Hence, the kinetic energy flux can be 
expressed by Eq. ( pOf ) in terms of the effective mode-to-mode transfer ^ UIl (k|p|q). The kinetic energy flux II^> (the 
energy transfer rate from the it-sphere to outside the same sphere) in MHD will also be given by the same equation, 



K<(K)= ]T E r u (k| P |q). (105) 

\p\<K \k\>K 

The magnetic energy flux U^(K) is defined as the rate of energy lost by the ^-sphere to the modes outside the 
6-sphere. Since ljl bb (k\p\q) is the mode-to-mode transfer from b(p) to b(k), Il^^K) can be obtained by summing 

% (k|p|q) 

over every mode p inside the &-sphere and every mode k outside the fe-sphere, i.e., 

n£(#)= E E rtlplq). (106) 
|p|<K |k|>if 



As in the kinetic energy flux, the circulating transfer Ya will not contribute to the magnetic energy flux in Eq. (106). 
We briefly explain this as follows: let us take the mode b(q) to be outside the 6-sphere. Then if the mode b(k) loses 
Ya to mode b(p) then it will also gain Ya from mode b(q), and hence Ya will not contribute to the magnetic energy 
flux. Similar arguments apply when b(q) is inside the 6-sphere. Therefore, we can write the magnetic energy flux in 
Eq. ( |i06|) as 

n£(*)= E E A k lp|q)- (107) 

|p|<K |k|>K 



where ^ &b (k|p|q) is given by Eq. @ of Section |VII B| . 



Similarly, the energy transfer between the kinetic energy and magnetic energy can be described by four types of 
fluxes. We shall define these fluxes below. 

There is a transfer of energy from a u-sphere of radius K to the 6-sphere of the same radius. The rate of loss of 
energy from the u-sphere to the corresponding 6-sphere can be calculated by by summing ^. bll (k|p|q) over every mode 
inside the fe-sphere and the u-sphere. Following the definition of the effective shell-to-shell transfer T^ n we can define 
the effective flux from the if-sphere to the 6-sphere as 

k<(k)= E E ^"( k ip|q)- ( 108 ) 

|p|<K|k|</f 
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by summing the effective mode-to-mode transfer, $ bu (k\p\q), from u(p) to b(k), over the modes inside the w-sphere 
and the ^-sphere. 

There is a transfer of energy from the u-sphere to the magnetic modes outside the corresponding ^-sphere. The rate 
of loss of the the energy from the w-sphere to the modes outside the 6-sphere can be obtained by summing ^?. b "(k|p|q) 
over every mode inside the it-sphere and every mode outside the 6-sphere. The corresponding effective flux can be 
calculated using 

K>(«)= E E ^"( k ip|q)- ( 109 ) 

|p|<X \k\>K 

Similarly, there is a transfer of magnetic energy from a &-sphere to velocity modes outside the w-sphere. The rate 
of loss of energy from the ^-sphere to the modes outside the u-sphere can be calculated by 

nii(K)= E E ^ b (k| P |q). (110) 

|p|<if |k|>K 

There is a transfer of energy from modes outside the it-sphere to the modes outside the &-sphere. The rate of loss 
of energy by the modes outside the M-sphere to the modes outside the &-sphere can be obtained by 

n$£(#")= E E ^"( k lp|q)- ( m ) 

|p|>if |k|>K 

The total effective flux is defined as the total energy (kinetic+magnetic) lost by the if -sphere to the modes outside, 
i.e., 

u tot (K) = n«< (K) + n b < (K) + n h "< (k) + n b u < (k). (in) 



A schematic illustration of the effective fluxes defined in Eqs. ( |105[ )-(111) was given in Fig. 10. If the circulating 
transfer is zero, then each of the effective fluxes will be the same as the corresponding total fluxes. In a subsequent 
paper we will present the results of the detailed study of these fluxes in numerical simulations of 2-D MHD turbulence. 



IX. CONCLUSION 

In literature we find a description of energy transfer rates from two modes in a triad to the third mode, i.e., from 
u(p) and u(q) to u(k). Here we have constructed new formulae to describe energy transfer rates between a pair of 
modes in the triad (mode-to-mode transfer), say from u(p) to u(k). The third mode in the triad acts as a mediator 
in the transfer process. 
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The mode-to-mode energy transfer in our formalism can be expressed as a combination of an "effective transfer" 
and a "circulating transfer" . In Appendices A-C we have shown by imposing galilean invariance, symmetry and some 
other constraints that the circulating transfer should be zero. A completely general proof is not available at present. 
However even if the circulating transfer happens to be nonzero, it will not result in a change of modal energy, since 
the amount of circulating transfer gained by mode k from mode p is also lost by mode k to mode q (see Figs. 3, 8, 
9). Only the effective transfer is responsible for modal energy change. As the circulating transfer does not have any 
observable effect on the energy of the modes, it may be correct to ignore it from the study of energy transfer. 

Using the notion of effective mode-to-mode transfer and the circulating transfer we defined effective shell-to-shcll 
energy transfer rates and energy fluxes in fluid turbulence [see Eq. (|27|) and Eq. J30|)1 and in MHD turbulence [see 



Eqs. (102 )- ( J104| ) and Eqs. (|1 5| ) -(111)]. Some of these energy transfers can be obtained only using the "mode-to-mode" 
energy transfer formulae. In a companion paper fl7| we present the results of our numerical study of the effective 
shell-to-shell energy transfer rates and energy fluxes in 2-D MHD turbulence. 
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APPENDIX A: DERIVATION OF THE CIRCULATING TRANSFER BETWEEN THE VELOCITY 

MODES IN A TRIAD 



In Section IV we showed that the mode-to-mode transfer is given by TjC 1 = $ u + Xa- We interpreted Xa as a 
circulating transfer and $ uu [Eq. (p"2])] as an effective mode-to-mode transfer. In this appendix we shall determine the 
circulating transfer Xa for fluid and MHD turbulence, using symmetry considerations. In the first section we shall 
discuss Xa in fluid turbulence and in the next section we will discuss the MHD turbulence. 

1. Xa in fluid turbulence 

The energy transfer between pair of modes should be invariant under a rotation (because it is a scalar) , and galilean 
transformation, and it should be finite for all values of k, p, q, u(k), u(p), u(q). We will determine the form of the 
circulating transfer by imposing these and a few other symmetry constraints. 

We construct a scalar expression for Xa dependent on k, p, q, u(k), u(p), u(q), and which is cubic in u and 
linear in k, p or q, to satisfy the for dimensional reasons. We write the general expression for such a scalar as 

X A = i?e([k.u(q)]{[u(k).u(p)]a + [u(k).u(q)]/3 + [u(p).u(q)] 7 + 

[u(k).u(k)]S + [u(p).u(p)]i/ + [u(q).u(q)M + 
[k.u(p)]{[u(k).u(p)] ?? + [u(k).u(q)]C + [u(p).u(q)]e + 

[u(k).u(k)]/s + [u(p).u(p)]A + [u(q).u(q)]0} + 
[p.u(k)]{[u(k).u(p)]^ + ... + [u(p).u(q)]u; + ... + [u(q).u(q)]x} + 
[p.u(q)]{[u(k).u(p)]r + ... + [u(p).u(q)]<£ + ... + [u(q).u(q)]p} + 
[q.u(k)]{[u(k).u(p)]a + ... + [u(p).u(q)]£ + ... + [u(q).u(q)]?} + 

[q.u(p)][[u(k).u(p)]7r + [u(k).u(q)]p+ ... + [u(q).u(q)]e}) (A.l) 
where all the coefficients a, (3, etc. are atmost non-dimensional functions of the wavenumbers k, p, q and of u(k), 



u(p), u(q). 'Re' in the expression in Eq. (A.l) indicates the real part of the expression. We should have infact 
written Xa as a linear combination of the real and the imaginary parts of the entire expression in the above equation. 
However, we will only consider the real part explicitly. The imaginary part can be treated in the same manner and 
the following conclusions, although explicitly stated for the real part will also be valid for the imaginary part. The 
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most general expression should also contain terms of the kind i?e[k.u(q)]i?e[u(k).u(p)], i?e[k.u(q)]7m[u(k).u(p)], 
and 7m[k.u(q)]im[u(k).u(p)]. However, we will explain below that such terms should not be present in AT a on 
account of galilean invariance. Hence, from the outset we do not include such terms. The task now is to determine 



the coefficients in Eq. (A.l) 



The interactions involves all three modes of the triad, i.e., if the Fourier coefficients of any of the wavenumbers is 
zero, interaction between the modes is turned off. The mode-to-mode transfer should respect this fundamental feature 
of the triad interactions. Therefore, if u(k), u(p) or u(q) approach zero, Xa should approach zero as well. We 



will argue below that a,/?, etc., in Eq. (A.l) are independent of the magnitude of u(k), u(p), u(q). By inspecting 



Eq. (A.i) for Xa, hence we find that a few terms in the expression are independent of the magnitude of atleast one 
of the Fourier components. For example, the term [k.u(q)] [u(k).u(q)]/3 is independent of u(p). Such terms should 
be dropped from the expression of circulating transfer. After dropping such terms, we obtain 

X A =#e([k.u(q)][u(k).u(p)]a+ [k.u(p)][[u(k).u(q)]C + 
[p.u(k)][[u(p).u(q)] W + [p.u(q)][[u(k).u(p)]r + 

[q.u(k)][[u(p).u(q)]f+[q.u(p)][[u(k).u(q)]p) (A.2) 



The requirement of finiteness of the mode-to-mode transfer imposes restrictions on the form of Xa- In Eq. (A.l) 
the dimensional terms are [k.u(q)][[u(k).u(p)], [k.u(p)][[u(k).u(q)], etc. and they are finite for all values of k, p, q, 
u(k), u(p), u(q). The coefficients a, (3, etc. are 'dimensionless scalars'. Therefore they can be written as 

= f\ k p k p k u (P) k u (p) u(k).u(p) u(k).u(p) 

a I[ kp> kq "^^'k\u(p)\'k\u{<^)\''" , |u(k)||u(p)|' |u(k)||u(q)|'- J 1 ' j 

where the arguments of the function have been non-dimensionalised. Consider triads in which one of the wavenumbers, 
say q, tends to zero. Then the terms like k.p/fcg with q in the denominator diverge. Similarly if one of the Fourier 
coefficients, say u(q), is equal to zero, then the arguments k\u(q)\ ano - |u(k)|'|"^(q)| 1 e ^ c - diverge. Hence terms which 
depend on the magnitude of vectors must not appear in the expression of the dimensionless scalars a, (3, etc. Only the 
arguments which depend on the angles between vectors, e.g., k\u(p)\ ' \u(k)\\u(p) \ ' snou ld t> e included. Therefore, all the 
dimensionless coefficients should be a function of only these arguments, i.e., a = /[^, k\u(p)\ > |u(k)|'|u(p)| ' "']' 
We now write the coefficient a, etc. in the form 

a = «(°) + + a (2) + a (3) + a (4) + a (5) (A 4) 
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where a' ' is a constant, and 



a(2) = ?) tuW a ( 2 ) Wq) Q (2) P^k) + a (2) P_u(q) Q ( 2 ) q-W + ^ f^P) (A>(J) 

fe|u(p)| fc|u(q)| p|u(k)| p|u(q)| g|u(k)| g|u(p)| 

Q (3) _ a (3) u(k)-u(p) (3 ) u(k).u(q) (3) u(p).u(q) 

" ai |u(k)||u(p)| +a2 |u(k)||u(q)| +a3 |u(p)||u(q)|- ^ 

The quantities a^ 2 \ a' 3 -* are 'linear' in the scalars of the type k.p/fcp, k.u(p) /fc|u(p)|, and u(k).u(p)/|u(k)||u(p)| 
respectively. The quantity contains terms like (k.p) 2 /k 2 p 2 , [k.u(p)] 2 /fc 2 ([u(p).u(p)] , etc. which are higher order 



in these scalars. All the other coefficients £,7, w, etc. in Eq. (A.f) can be similarly divided into linear and non-linear 
terms. 



We now study the consequences of galilean invariance of on Eq. ( A.l ). Consider two frames of reference S and 
S' with S moving with respect to S' with velocity Uo- We denote the Fourier coefficient of velocity in frame 5* and 
S' by u(k) and u'(k) respectively. The relationship between u(k) and u'(k) is given by 

u'(k) = u(k) exp 4(k ' Uo)t +U (5(k). (A.8) 

The second term is non-zero only for the k = modes. In the subsequent discussion relating to galilean invariance we 
can drop this term for sake of simplicity without affecting the derivation. The transformation rule for scalars of the 
type k.p, k.u(p), u(k).u(p) between the two frames of references S and S' can be deduced from the transformation 
property of u(k). 

k'.u'(q') = k.u(q) exp l(q Uo)t (A.9) 

u'(k').u'(p') = u(k).u(p) exp ,; ( k + p >- Uo * (A.10) 

It follows from the above two transformation rules that k.u(q) and u(k).u(p) are not invariant under a change of 
reference frame. However, a product of these quantities can be seen to be invariant 

[k'.u'(q')][u'(k').u'(p')] = [k.u(q)][u(k).u(p)] exp l(k+p+q) - Uot 

= [k.u(q)][u(k).u(p)] (A.ll) 

using k+p+q=0. However, 
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i?e[k'.u'(q')]i?e[u , (k').u'(p')] + i?e[k.u(q)]i?e[u(k).u(p) 



(A.12) 



and 



ite[k\u'(q')]Iro[u'(k').u / (p')] ^ i?e[k.u(q)]/m[u(k).u(p)], 



/m[k'.u'(q')]/m[u'(k').u'(p')] + iro[k.u(q)] Jm[u(k).u(p)] 



(A.13) 



(A.14) 



Thus, we have shown that the terms of the kind given in the above three equations are not galilean invariant. 

Clearly, the invariance under Galilean transformation demands that all of the vectors u(k), u(p), u(q) appear 
either once, twice, or n times in the scalar. A scalar of the type [k.u(q)][u(k).u(q)] having two or more repeated 
wavenumber indices is not invariant since the term appearing in the exponential is not zero for such a scalar. The 
transformation properties of any quantity formed from the scalars like k.p, k.u(p), u(k).u(p) can be obtained 
in a similar manner. If the quantity is not galilean invariant then it should not be included in the expression for 
mode-to-mode transfer. 



We have shown above that scalars of the type k.u(q) and u(k).u(p) are not Galilean invariant [Eqs. (A. 9) and 
( |A.10| )1. It follows that ,/?( 2 \ etc., a , etc. are not Galilean invariant. To preserve the Galilean invariance 



of Xa, the coefficients a*- 2 -* and a^ 3 ) should be dropped from Eq. (A4). However, the nonlinear quantity a*- 4 -* can be 
constructed in a Galilean invariant form. Currently, we have dropped these higher order terms. Our derivation of Xa 
therefore has the limitation of being strictly valid only upto the linear order. A complete proof containing all order 
of terms in Aa is beyond the scope of this thesis and is left for future work. After dropping the terms , a*- 3 -* , , 



the Eq. (A.4) reduces to 



a = a<°> 



(A.15) 



Using equations ( A.4 ) and ( A.15| ) we now write X in Eq. (A2) as 

X A = Jte([k.u(q)] 
[k.u(p)] 
[p.u(k)] 

[P-u(q)] 
[q.u(k)] 

[q-u(p)] 



u(k).u(p)]{ ao + «i^+« 2 -^ + «3 — } 

u ( k).u(q)] {C o + C 1 ^ + C 2 ^ + C3^} + 

M Mlf k.p k.q p.q 

u(p).u(q)]{^ + +^ 2 — +w 3 — } 

, n , m k.p k.q p.q 

u(k).u(p) {t + ti— +t 2 - hr 3 } + 

kp kq pq 

u(p).u(q)]U + ^ + & — } + 

kp kq pq 

u(k).u(q)]{p + Pl — + p 2 — +p3 — }) 



(A.16) 
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We have dropped the superscript '1' from , a^- Using the incompressibility relationship q.u(q) = we 
obtain p.u(q) = — (k + q).u(q) = — k.u(q). Therefore, fourth term in the above equation can be absorbed into the 
first term. Similarly, using the relationships q.u(p) = — k.u(p) and q.u(k) = — p.u(k), the last term can be combined 
with the second term and the fifth term combined with the third term. Now, we can write 

X A = i?e([k.u(q)][u(k).u(p)]{a + + a 2 ^ + a 3 — } + 

kp kq pq 

[k.u(p)][u(k).u(q)]{Co + Cl^ + C2^ +C3 — } + 

kp kq pq 

tp.u(k)][u(p).u(q)]{wo + + + lo^}) (A.17) 

by making the replacements (a — r) — > a, (( — p) —> C an d (oj — £) — > 

The quantity X A represents the energy circulating in the anti-clockwise direction, i.e., X\ flows along p — » k > 



q — > p [see Fig. 3 in Section IVB]. We will now denote the circulating transfer along p — > k by A(k|p|q), transfer 
along k — > q by A(q|k|p) and so on for the circulating transfer between each pair of modes. The above X A is the 
circulating transfer from p to k, i.e., X(k|p|q|) = Then by symmetry the circulating transfer from q to k can 
be written as 

X(k|q|p) = i?e([k.u(p)][u(k).u(q)]{a + £*i-^ + a 2 ^- + a 3 — } + 

f£Q fvj) J)Q 

[k.u(q)][u(k).u(p)]{Co + Ci^ + C2^ + Cs^} + 

[q.u(k)][u(p).u(q)m + U} ^+ + ^3^}), (A.18) 

and the circulating transfer from p to k is 

A(p|k|q) = fle([p.u(q)][u(k).u(p)]{a + + a 2 — + a 3 ^} + 

kp pq kq 

[p.u(k)] [u(p) .u(q)] {Co + Ci ^ + (2 — + Cs ^ } + 

kp pq kq 

[k.u(p)][u(k).u(q)]{wo + + ^2— + ^3^}) (A.19) 

From Fig. 3, AT(k|p|q) flows along the anti-clockwise direction while X(k|q|p) and A(p|k|q) flows along the clockwise 
direction. The former should therefore be equal and opposite to the latter two transfers, i.e., 

X(k|p|q)+X(k|q|p)=0, (A.20) 

and 
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X(k|p|q)+A(p|k|q) =0. 



(A.21) 



Substitution of expressions A(k|p|q), X(k|q|p), A(p|k|q) [Eqs. ( |A.17| )-( [09| )1 into the first of these two relationships 
yields 

[k.u(q)][u(k).u(p)]{(a + Co) + (a x + (2)^ + (a 2 + Ci)^ + (a 3 + Ca)^} 
[k.u(p)][u(k).u(q)]{(Co + ao) + (Ci + « 2 )^ + (Ca + + (Ca + «a)^} 

[p.u(k)][u(p).u(q)]{(w! -w 2 )^- + {uj 2 - wi)^} = 0, 

(A.22) 



and 



[k.u(q)][u(k).u(p)]{(a 2 - a 3 )^ + (a 3 - a 2 ) — } 
-[k.u(p)][u(k).u(q)]{(Co + wo) + (Ci + (C2 + ^)^ + (Ca + ^2)^} 

-[p.u(k)][u(p).u(q)]{(a;o + Co) + (wi + Ci)^ + (u* + Ca)^ + (w 3 + C 2 )^}= 0. 



(A.23) 



Each of the three terms in the two equations above has a different functional dependence on u(k), u(p), u(q). If 
we only rotate u(k) about k without changing its magnitude then each of the terms will change by a different factor 
depending on the direction of u(p), u(q) which can be independently chosen. Since the sum of three terms should 
be zero and each term can change a different factor, each of the terms should be individually zero. In addition, 
this should hold for all wavenumber triads. Therefore each term inside the curly bracket should also be individually 
zero. This condition gives the the following relationships between the constants: ctQ + Co = 0, a± + C2 = 0, a 2 + Ci = 



0, a 3 + C3 = 0, uji - LU2 = from Eq. ( A.22 ), and uj + Co = 0, a 2 — a 3 = 0, £i + u>i = 0, C2 + w 3 = 0, C3 + ^2 = from 



Eq. ( |A.23| ). 

We now consider the following geometry - the wavenumber triad form an isosceles triangle with |p| = |q|, and 
the Fourier modes u(p) and u(q) are perpendicular to the plane of k, p, q and are equal in magnitude. For this 
geometry the clockwise and the anti-clockwise direction are equivalent. There is no preferred direction along which 
the circulating transfer may flow. Therefore, for the geometry shown the circulating transfer should be zero. We may 



now write the circulating transfer [Eq. (A. 17)] for this geometry as 



X = ite([p.u(k)][u(p).u(q)K(«i + u 2 )^ + W3 — + a}) 

Kp pq 
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= (A.24) 

With u(p) and u(q) fixed perpendicular to the plane of the triad, X should be zero for all wavenumber triads which 
form an isosceles triangle. For X to be zero for every such geometry, each term inside the curly bracket should be 
equated to zero giving the additional relationships: lo\ — — cj2,^3 = 0,c<jo = 0. 

We have obtained 17 relationships between the constants a, (, etc. These equations may be solved to get the 
values of all these constants. It can be verified that the only solution of these equations is a's= 0, C' s — 0, oj's= 0, 
a = b = c = 0. 

Therefore we get X = 0. It follows that the mode-to-mode transfer T/l uu = fl uu . 

We have thus determined the mode-to-mode transfer. The proof used rotational invariance, galilean invariance, 
finiteness of mode-to-mode transfer, and symmetry of X with respect to k, p, q, u(k), u(p), u(q) . However, in the 
series we have taken only linear terms in k.p, k.u(p), u(k).u(p). A complete proof containing all orders of terms is 
beyond the scope of this thesis. It may be possible to obtain a completely general proof by finding some symmetries 
which we have ignored in our proof. Otherwise a novel approach may be required. 



2. X A in MHD turbulence 

The circulating transfer between the velocity modes in MHD can also be shown to be equal to zero to linear 
order terms in k.p, k.u(p), u(k).u(p) . The arguments leading to this are nearly the same as that for the fluid case. 
However, we additionally need to show that the general expression for the circulating transfer X& for MHD should 



not contain b(k), b(p), b(q) and can also be given by Eq. (AT) 

In the previous section we showed that each term in A a should depend on u(k), u(p), u(q) — this followed from 
the fundamentally three-mode feature of the interactions. The expression for the combined energy transfer to u(k) 
from u(p) and u(q) is identical for the fluid and for the MHD case. Hence, in the case of MHD also, each term in 
the circulating transfer between the velocity modes should be dependent on u(k), u(p), and u(q). Thus, the general 
expression for the circulating transfer between the velocity modes in MHD is the same as that for Fluid turbulence. 



The rest of the arguments leading to Xa = in MHD are identical to those in Section A 1 
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APPENDIX B: DERIVATION OF CIRCULATING TRANSFER BETWEEN THE MAGNETIC MODES 

IN A TRIAD 



In Section VII B we had shown that the mode-to-mode transfer from b(p) to b(k) can be written as T/t bb = $ bb + Y&, 



where $ bb given by Eq. ( Jl^ ) is the effective transfer between the modes and Ya is the circulating transfer. In this 
section we will derive the expression for the circulating transfer Ya- The derivation for Ya is similar to the derivation 
of Xa for fluid turbulence given in Appendix A. 

Since the mode-to-mode transfer is a scalar, the circulating transfer Ya will also be a scalar with a dependence on 
the wavenumbers k, p, q and the Fourier coefficients u(k), u(p), u(q), b(k), b(p), b(q). It should be linear in k, 
p or q and cubic in the Fourier coefficients u(k), u(p), u(q), b(k), b(p), b(q) for dimensional reasons. 

The most general form of Ya that satisfies these properties can be written as 

Y A =Yi 1) +Y^+Yi 3) +Yi 4) (B.I) 

where 

Y™ = J Re([k.u(q)]{[u(k).b(k)]... + [u(k).b(p)]... + [u(k).b(q)] ... + 

[u(p).b(k)]... + [u(p).b(p)]... + [u(p).b(q)]... + 

{[u(q).b(k)]... + [u(q).b(p)]... + [u(q).b(q)]...} + 
[k.u(p)]{[u(k).b(k)]... + .... + [u(k).b(q)]... + .... + [u(q).b(k)]... + .... + [u(q).b(q)] ...} + 
[p.u(k)]{[u(k).b(k)]... + .... + [u(p).b(q)]... + .... + [u(q).b(p)]... + .... + [u(q).b(q)] ...}) 

(B.2) 

Yf = fle([k.u(q)]{[b(k).b(k)] + [b(k).b(p)]a + [b(k).b(q)]... + 
[b(p).b(p)] + [b(p).b(q)]... + [b(q).b(q)]...} + 
[k.u(p)]{[b(k).b(k)]... + .... + [b(k).b(q)]/3 + .... + [b(q).b(q)]...} + 
[p.u(k)]{[b(k).b(k)]... + .... + [b(p).b(q)] 7 + .... + [b(q).b(q)]...» 

(B.3) 



Y^ = i?e([k.b(q)]{[u(k).b(k)]... + [u(k).b(p)]* + [u(k).b(q)].. 
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[u(p).b(k)]Tj + [u(p).b(p)]... + [u(p).b(q)]... + 

{[u(q).b(k)]... + [u(q).b(p)]... + [u(q).b(q)]...} + 
[k.b(p)]{[u(k).b(k)]... + .... + [u(k).b(q)] M + .... + [u(q).b(k)]^ + .... + [u(q).b(q)]...} 
[p.b(k)]{[u(k).b(k)]... + .... + [u(p).b(q)ty + .... + [u(q).b(p)] X + [u(q).b(q)] ...}) 



(B.4) 



Y A 4) = i?e([k.b(q)]{[u(k).u(k)]... + [u(k).u(p)]0 + [u(k).u(q)] ... + 
[u(p).u(p)]... + [u(p).u(q)]... + [u(q).u(q)]...} + 
{[u(q).b(k)]... + [u(q).b(p)]... + [u(q).b(q)]...} + 
[k.b(p)]{[u(k).u(k)]... + [u(k).u(q)]0 + .... + [u(q).u(q)]...} + 
[p.b(k)]{[u(k).u(k)]... + [u(p).u(q)]e + .... + [u(q).u(q)]...]) 



(B.5) 



Following Appendix A we consider explicitly only the real part, which is denoted by Re in the above equations. 
If we take a linear combination of the real and the imaginary parts instead, the derivation and the conclusions 
regarding Ya will not change. In addition, terms of the i?e[k.u(q)]i?e[u(k).u(p)], i?e[k.u(q)]i?e[u(k).b(p)], etc. are 
not included as they violate galilean invariance - we had explained this in Appendix A in context of The three 
dots in the above equations in front of some of the terms indicate that a coefficient multiplies each of these 
terms but which we have not shown explicitly as they do not enter most of the future discussion. The coefficients 
a, (3, etc. in the above expressions for Y^\.. Y^ are non-dimensional functions of the wavenumbers k, p, q and 
of u(k), u(p), u(q), b(k, b(p), b(q). To determine Ya we should now determine these coefficients. Note that 



besides the terms of the kind [k.u(q)]{...}, [k.u(p)]{...}, [p.u(k)]{...} in Eqs. (B.2)-( B.5 ), there are also terms of the 
kind [p.u(q)]{...}, [q.u(p)]{...}, [q.u(k)]{...}. However, using the incompressibility constraints k.u(k) = and the 
constraint k + p + q = 0, we have absorbed the latter three terms into the former ones. Now we will impose a few 
restrictions on the form of Ya • 

The interactions between the magnetic modes fundamentally involve three-modes — the combined energy 
transfer to a magnetic mode, say b(k), depends on the Fourier coefficients at all three wavenumbers (k, 
p, q). For example, the combined energy transfer to b(k) from b(p) and b(q), given by S ,hb (k|p, q) = 
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-i?e(i[k.u(q)][b(k).b(p)] + z[k.u(p)][b(k).b(q)]), depends on u(p), u(q), b(p), b(q), and b(k). The mode-to- 
mode transfer should not break the fundamentally three- mode feature of these interactions, i.e., the mode-to-mode 
transfer should have Fourier coefficients from all wavenumbers k, p, q. Hence we will drop all those terms from 
to Y^ 1 which are independent of atleast one of the wavenumbers k, p, q. Then, we get 

y« = J?e([k.u(q)] {[u(k).b(p)]... + [u(p).b(k)]...} +[k.u(p)] {[u(k).b(q)j... + [u(q).b(k)] ...} + 

[p.u(k)] {[u(p).b(q)]... + [u(q).b(p)]...» (B.6) 

Y A 2) = i?e([k.u(q)][b(k).b(p)]«+ [k.u(p)][b(k).b(q)]/3+ [p.u(k)][b(p).b(q)] 7 ) 

(B.7) 

y| 3) = J?e([k.b(q)]{[u(k).b(p)]<5 + [u(p).b(k)]r 7 } + 
[k.b(p)]{[u(k).b(q)] M + [u(q).b(k)]i/} + 

[p.b(k)]{[u(p).b(q)]^ + [u(q).b(p)] X }) (B.8) 
y| 4) = fle([k.b(q)][u(k).u(p)]0+ [k.b(p)][u(k).u(q)]0+ [p.b(k)][u(k).u(k)]...) 

(B.9) 

We put the following additional restriction on the form of Ya- If none of the magnetic modes in the triad are 
gaining/losing energy to the other two modes, i.e., if the combined energy transfer to each one of them is equal to 
zero, then even the circulating transfer should be equal to zero. We observe the following additional feature of the 
interactions. The combined energy transfer to every magnetic mode in the triad from the other two magnetic modes 
is equal to zero, if any two of b(k), b(p), b(q) are zero (see Eq. (|36|). The mode-to-mode transfer should also respect 
this feature of the interactions, i.e., the circulating transfer should be zero if {b(k), b(p)=0}, or {b(k), b(q)=0}, 
or {b(p), b(q)=0}. Now, let us put b(k)=0, b(p)=0 in Y^ to Y^\ It can easily shown that Ya will be zero, for 
arbitrary values u's and b's, only if the all coefficients in Y^' and Y^' are themselves equal to zero. Hence we drop 



Y^ 1} and Yj 4) from Y A [Eq. (fy])]. We get 



and 
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Y A - Y| 2) + Y^ 3) (B.10) 



r| 2) = i?e([k.u(q)][b(k).b(p)]a+ [k.u(p)][b(k).b(q)]/3 + 

[p.u(k)][b(p).b(q)] 7 ) (B.11) 

r| 3) = ite([k.b(q)]{[u(k).b(p)]J + [u(p).b(k)M + 
[k.b(p)]{[u(k).b(q)]/i + [u(q).b(k)]i/} + 

[p.b(k)]{[u(p).b(q)]V + [u(q).b(p)] X }) (B.12) 

In Appendix A we had considered the derivation of circulating transfer between the velocity modes. By imposing 
the condition of finiteness of the circulating transfer we had concluded that the coefficients should be independent 
of the magnitude of the wavenumbers and the Fourier coefficients u(k), etc. On imposing the constraint that Y\ 
should be finite and following the reasoning in Appendix A, we can conclude that the coefficients a, (3, etc. should 
be independent of the magnitude of the wavenumbers and of the Fourier coefficients of the velocity and the magnetic 

fiplrk TWpfnrp wp mav writp thp rnpffiripnt rv — fT kp ku (p) kb (p) "(k).u(p) b ( k )- b (p) "(k).b(p) 1 

tieias. inereiore, we may write trie coefficient a - j^^, j\b{pj\^ -> \u(k)\\u( P )\ > |fc(fc)||f>( P )| ' -' \u(k)\\b(p)\ > -J- 

The other coefficients can also be written in the same manner. 
We will now write the coefficients a, (3, etc. as 

a = a<°> + + a (2) + a (3) + a (4 \ (B.13) 



where is a constant, and 



( 2 ) (2) k.u(p) (2) k.u(q) ( 2 ) P-u(k) (2) p.u(q) (2) q.u(k) (2) q.u(p) 

fc|u(p)| fc|u(q)| 3 p|u(k)| p|u(q)| 5 g|u(k)| 6 g|u(p)| 
(2) kb(p2 (2) q.b(p) 
" 7 fc|b(p)| + - + ai2 g |b(p)|' (R15) 



(3) = (3) u(k).u(p) (3) u(k).u(q) (3) u(p).u(q) 

1 |u(k)||u(p)r 2 |u(k)||u(q)r 3 |u(p)Hu(q)l 
(3) b(k).b(p) (3) b(k).b(q) (3) b(p).b(q) 

"4 i. „ mi, / m T « 5 T « 6 



|b(k)||b(p)| 5 |b(k)||b(q)| b |b(p)||b(q)| 
(3) u(k).b(k) (3) u(k).b(p) (3) u(q).b(p) 

" 7 |u(k)||b(k)| +a8 |u(k)||b(p)| + - + ai5 |u(q)||b(p)r ^ 

where a^a^af , etc. are all unknown constants. The quantities c^ 1 ), a^ 2 \ are 'linear' in the 
scalars of the type k.p/fcp, k.u(p)/fc|u(p)|, k.b(p)/fc|b(p)|, u(k).u(p)/|u(k)||u(p)|, b(k).b(p)/|b(k)||b(p)|, and 
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u(k).b(p)/|u(k)||b(p)|. The quantity contains higher order terms in these scalars like (k.p) 2 /fc 2 p 2 , 

[k.u(p)] 2 /fc 2 ([u(p).u(p)], etc. All the other coefficients £,7,0;, etc. in Eqs. (B.11)-(B.12) can be similarly divided 
into linear and non-linear terms. 

The mode-to-mode transfer should be invariant under a galilean transformation. In Appendix A, we had given 
the transformation rules for the scalars of the type k.p, k.u(p), u(k).u(p) and for a combination of these like 
[k.u(p)] [u(k).u(q)] . For the present case we also need to consider the transformation properties of the magnetic 
field. MHD equations hold for low velocities and only quantities upto the order of v/c are included in the equations. 
Under this approximation the magnetic fields are equal in the two frames S and S' with a relative velocity Uo. Then, 
the Fourier coefficient b(k) in the frame S and b'(k) in the frame S' are related as 



b'(k) =b(k)exp i ( ku °)* 



(B.17) 



This differs from the galilean transformation of u(k) by the term Uo<5(k) only. Using the above relationship between 
between b(k) and b'(k) the transformation rules for k.b(q), b(k).b(p), u(k).b(p) are 



k'.b'(q') = k.b(q) ex P 4(q U °)*, 



(B.18) 



b'(k').b'(p') = b(k).b(p) exp^+P}-" *, 



(B.19) 



u'(k').b'(p') = u(k).b(p) exp l(k+p)u °* . 



(B.20) 



These scalars are not galilean invariant (in Appendix A, we had likewise shown that the scalars of the type k.u(p), 
u(k).u(p) are not galilean invariant). It follows that a^ 2 \f3^ 2 \ etc., a^ 3 \ (3^- 3 \ etc. are not invariant under a galilean 
transformation. To preserve the galilean invariance of Ya, the coefficients a^ 2 \ cS 3 \ (3^ 2 \ f3^ 3 \ and other such 



coefficients should be dropped from Eq. (B.13). The nonlinear quantity f3^\ etc. can be constructed in a 

galilean invariant form. Even then we will drop this term. We had also dropped this term from the expression for 
X/\ in Appendix A. Our derivation for Y"a, like the derivation of A a in Appendix A, is is thus valid only upto the 



linear orders. After dropping the coefficients offl, a^ 3 \ and other coefficients of this type, Eq. (B.13) reduces to 



a = a l°>+« (1 > 



(B.21) 



Using Eqs. ( |B.2lD and flB.14| ) we can now write and in Eqs. Q.ll| )-( |B~1^ ) as 
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ll 2) = i?e([k.u(q)] [b(k).b(p)]{«W + ai % + a 2 ^ + a 3 ^} + 
[k.u(p)][b(k).b(q)]{/3(°) +0^+^ + folg} + 
[p.u(k)][b(p).b(q)]{ 7 (°)+ 7l ^ +72 k ^+ 73 ^} ) (B.22) 



Ret fk bfallfufkl bCnll-f^ ) + <^ ^£ 4 






[k.b(q)][u(p).b(k)]{ ?7 ( )+ f?1 ^4 






[k.b(p)][u(k).b(q)]{ M (°) +Ml k ^4 


kq 


+ ^} ■+ 


[k.b(p)][u(q).b(k)]{^) + ^^4 


-2^ 


+ -3^} + 


[p.b(k)][u(p).b(q)]{^ )+^^4 






[p.b(k)][u(q).b(p)]{ X (°)4Xiif 


^ 2 ^ ^3 P p q } 



(B.23) 



where the superscript '1' has been dropped from a' 1 -', a' 2 \ etc. 



The quantity Ya = represents the energy circulating in the anti-clockwise direction, i.e. Ya flows along 



p — > k — > q — > p [see Fig. 8 in Section VII B]. We will now denote the circulating transfer along p — » k by I^(k|p|q), 



transfer along k — » q by Y/^(q|k|p) and so on for the circulating transfer between each pair of modes. We take 



Y A = Y A > + Yp from Eqs. (|B.22|) - (|B.23|) to be the circulating transfer from p to k, i.e., Y A (k|p|q|) = Fa- In the 



same manner >A(k|q|p)) = (k|q|p) + Y^ (k|q|p) can be written as 



and 



>l 3) (k|q|p) 



i?e([k.u(p)][b(k).b(q)]{a( '4 




k.p 

«2^f H 


"«3^} + 


[k.u(q)][b(k).b(p)]{/3(°)4 








[q.u(k)][b(p).b(q)]{ 7 (°) 




~ '2 fep 


+ 73^} ) 


Re( [k.b(p)][u(k).b(q)]{5(°)4 








[k.b(p)][u(q).b(k)]{^(°)4 






+ 


[k.b(q)][u(k).b(p)]{ M (°)4 


Mi if + 


M2^H 


-M3^} + 


[k.b(q)][u(p).b(k)]{^°)4 




-2^4 




[q.b(k)][u(q).b(p)]{<£<°> + 






+ 


[qb(k)][u(p).b(q)]{ X (°) 


+ Xi^ 


i k.p 

+ X2Tf 


+ X3^} ) 



(B.24) 



(B.25) 
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and the circulating transfer from p to k, Y" A (p|k|q)) = F A (p|k|q) + Y" A (p|k|q) can be written as 

Y A 2) (p|k|q) = Jfe(tp.u(q)][b(k).b(p)]{a(°) + + a 2 ^ + a 3 ±f} + 

[p.u(k)][b(p).b(q)]{/3(°) + (h$* + /3 2 ^ +/3 3 ^} + 
[k.u(p)][b(k).b(q)]{ 7 (°)+ 7l k ^+ 72 ^ +7 3^} ) (B.26) 



r A 3) (p|k|q) = Re( [p.b(q)][u(p).b(k)]{«5( ) 4 




r^^ 3 + 




[p.b(q)][u(k).b(p)]{^°)4 




[p.b(k)][u(p).b(q)]{ M (°)4 




"■^ pq 




[p.b(k)][u(q).b(p)]{^4 


kp 


PQ 


-3^} + 


[k.b(p)][u(k).b(q)]{^°) + 




rz pq 


+ 


[k.b(p)][u(q).b(k)]{x(°) 




X- 2 pq 


+ X3^} ) 



(B.27) 

YA(k|p|q) flows along the anti-clockwise direction while >A(k|q|p) and Y"A(p|k|q) flow along the clockwise direction. 
The former should therefore be equal and opposite to the latter two transfers, i.e., 

r A (k|p|q) + r A (k|q|p) = 0, (B.28) 

and 

y A (k|p|q) + r A (p|k|q) = 0. (B.29) 



Replacing the expression for F A (k|p|q) from Eqs. ( |B.22| )-( |R23| ) and for F A (k|q|p) from Eqs. ( |B.24| )-( [R25| ) into Eq. 
(B.28) we get 

[k.u(q)][b(k).b(p)]{(a(°) + p<®) + (ai + /? 2 )^ + (a 2 + Pi)^ + (a 3 + p 3 ^} + 

rZJ) rZQ J)Q 

[k.u(p)][b(k).b(q)]{(/3( ) + a (0) ) + (Pi + a 2 )^ + (p 2 + ai ^+ p 3 + a^} + 

[p.u(k)][b(p).b(q)]{( 7 i - 72 )^ + (72 " 7i)^ + 

kp kq 

[k.b(q)][u(k).b(p)]{(<$<°> + i/<°>) + (S 1 + M2 )^-P + (J 2 + + {Ss + M3 )E^} + 

/cp kq pq 

[k.b(q)][u(p).b(k)]{(r/°) + + i*)^H + (r /2 + ^i)^ + fa, + ^ 3 )^} + 

fcp /eg 

[k.b(p)][u(k).b(q)]{( M (0) + <5 (0) ) + (Mi + ^2)^ + (M2 + + ( M3 + <5 3 )i^} + 
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[k.b( P )][u(q).b(k)]{(^°) + + {v x + m )±v + { V2 + m)^ + + m)— } + 

kp kq pq 

[p.b(k)][u(p).b(q)]{(0(°) - X (0) ) + (01 ~ X 2 )^ + (02 - Xi)^ + (03 - X 3 ) — } + 

fep fcg pq 

[p.b(k)][u(q).b(p)]{(x (0) - (O) ) + (xi - 02)^ + (X2 - 0l)^ + (X3 - 03)^} = 



(B.30) 



and replacing the expression for lA(k|p|q) an d 3^(p|k|q) from Eqs. ( B.26| )-( B.27 ) into Eq.(B.29) we get 



[k.u(q)][b(k).b(p)]{(a 2 - a 3 )^ + (a 3 ~ a 2 ) — } + 

kq pq 

[k.u(p)][b(k).b(q)]{(/3(°) + 7 (0) ) + (01 + 72)^ + (02 + 73^ + 03 + 72 — } + 

/cp kq pq 

[p.u(k)][b(p).b(q)]{( 7 ( ) + 0<®) + ( 7l +0l)^ + (72 + + (73 + + 

[k.b(q)][u(k).b(p)]{(<^ - ^(0)) + (5i - + (S 2 - 773)^ + (<5 3 - ^2)^} + 

[k.b(q)][u(p).b(k)]{(r/°) - J(°>) + (r?a - fc)^ + fa, - f 3 )^ + ( m - <5 2 )^} + 

[k.b(p)][u(k).b(q)]{(p( > + 0(°)) + (mi + 00^ + (Ma + 3 )^ + (Ms + 2 ) — } + 

kp kq pq 

[k.b(p)][u(q).b(k)]{(i/°) + + ( Vl + Xi)^ + (vi + X3)^ + (V3 + X2) — } + 

kp kq pq 

[ P .b(k)][u(p).b(q)]{(0(°) + p(°>) + (0! + + (0 2 + p 3 )^ + (03 + M2)^} + 

[p.b(k)][u(q).b(p)]{( X (°) + V W) + (Xi + fl)^ + (X2 + ^3)^ + (X3 + ^2)^} = 



(B.31) 



Since each term in the above equation has a different functional dependence on u(k), u(p), u(q), b(k), b(p), b(q), 
each term should be individually zero. It follows that the constants in the above equations should satisfy the following 
relationships : from Eq. (B.30) we get + 0^ = O, ai +0 2 = 0, a 2 + X = 0, a 3 + 3 = 0, 71 - 72 = 0S {a) + /i (0) = 
0, Sx + it3 = 0, 6 2 + Mi = 0, <S 3 + m 3 = 0, 7?(°) + = 0, Vl + v 2 = 0, i l2 + v x = 0, r/ 3 + ^ 3 = 0, - x (0) = 0, 0i - X2 = 
0,02 -Xi = O,03-X3 = 0, and from Eq. (B.31) we get a 2 -a 3 = 0,p(°)+jW = 0,/3 1 + 7i = 0,/3 2 +73 = 0,/3 3 +72 = 
0, <y(°) - 77(0) = 0, <Ji - r?i = 0, <5 2 - V3 = 0, f a - m = 0, p (0) + (o) = 0, Mi + 0i = 0, M2 + 03 = 0, Ma + 02 = 0, + x (0) = 

0, ^2 + X3 = 0, M3 + 02 = 0. 

The above relationships can be arranged into the following groups: {a^ = —0^°' = 7^}, {5^ = —/i 1 - ' 1 = 0(°) = 
V {0) = -^ (0) = X (0) }> {"1 = -02 = 73}, {"2 = -01 = 7l = 72 = "3 = {#1 = -M2 = 03 = X3 = -vi = Vl}, 

{5 2 = -Mi = 0i = X2 = -va = V3}, = -M3 = 02 = Xi = ~ v x = m}- 

Let us now consider a triad with |p| = |qj, k.p = k.q and with the Fourier modes u(p), u(q), b(p), b(q) perpen- 
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dicular to the plane formed by this triad. We argued in Appendix A that the circulating transfer should vanish for 



such a geometry. For this geometry we can write from Eq. (B.22) and ( B.23) as 



Y A = iM[p.u(k)][b(p).b(q)]{ 7 (°) + (71 +72 )^ + 7 3^} + 



[p.b(k)] [u(p).b(q)]{<M°> + {4>l + 4>2)^§ + h^} + 



p<i 



[p.b(k)][u(q).b(p)]{ X ^ + (x 



X 2 ) kp + ^ 3 P pq} ) 



(B.32) 



The above equality should hold for all triads with |p| = |q| and k.p=k.q. Therefore each term inside the curly brackets 
should be individually zero. This leads us to the following additional relationships between the various constants in 
the above equation: T <°) = 0, 7l + 72 = 0, 73 = 0, </><® = 0, </> 2 + </> 3 = 0,^ = 0, x (0) = 0, X i + Xs = 0,X2 = 0, 
Similarly, by choosing |k| = |p|, and u(k), u(p), b(k), b(p) perpendicular to the triad plane we get the relationships 
a(°) = 0, ai = 0, a 2 + a 3 = 0, 6^ = 0,5i = 0, 6 2 + 5 3 = 0, r?(°) = 0, 771 = 0, 772 + m = 0. and by choosing |k| = |q|, and 
u(k), u(q), b(k), b(q) perpendicular to the triad plane we get f3^ = 0,@2 — 0,Pi + P3 = 0, //°) = 0,^2 = 0,/ii+/^3 = 
0, = 0, v 2 = 0, v x + y 3 = 0. 

Using these relationships the above groups take the following values: { a(°) — — /?(°) — 7 (°) = 0}, {5^ = 
- M (°) = = 77(0) = -1/(0) = x (°) = 0}, {ai = = 73 = 0}, {a 2 = -fa = 71 = 72 = a 3 = -fa = 0}, 

{Si = -fj, 2 = <p3 = X3 = -Vl = Vl = 0}, {^2 = "Ml = (f>l = X2 = -I>3 = V3 = S3 = A*3 = ~<h = ~Xl = v l = ~»?2}- 

Hence, some of the constants are equal to zero, and all the non-zero constants are related to each other by just one 
undetermined constant which we will denote by 'C, i.e., {52 = — Ml = ••• = C} 

Putting the values of these constants in the expressions for Y^ 2) and Y^ 3) in Eqs. (|B.22| )- (|B.23| ) we get 



Yl 2) = 



(B.33) 



Yl 3) = i?e([k.b(q)][u(k).b(p)]{C^ - C^} + [k.b(q)][u(p).b(k)]{-C^ + C^} + 
[k.b(p)][u(k).b(q)]{-C^ + C^} + [k.b(p)][u(q).b(k)]{C^ - C^} + 
[p.b(k)][u(p).b(q)]{C^ - C^} + [p.b(k)][u(q).b(p)]{-C^ + C^} ) 



(B.34) 



with a single undetermined constant. 
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For the determination of Xa in Appendix A, the above considerations were sufficient to determine However, 
in the present case we are still left with one undetermined constant. We will now impose the constraint of locality 
which we state as follows: In the limit k, q — * oo >> p, Ya should tend to zero. If remains finite then the noise, 
which is always inevitably present in any flow at the smallest scales, would interact directly with the large scales. 
Moreover, the locality constraint is also necessary in order that the circulating transfer is consistent with the effective 
transfer between b(k) and b(q) which goes to zero for in the above limit. Let us take u(k), u(q) to be perpendicular 
to the plane of the triad. Then we can write Ya as 



Y A 3) = J?e([k.b(p)][u(k).b(q)]{-C^ + C^} + [k.b(p)][u(q).b(k)]{C^ - C^}) 



kp 



PI 



which may be written in terms of the angles of the triangle formed by k, p, q [Fig. B.l] as 



Y A 3) =i?e( fc.b(p) [u(k).b(q)]{-C7ccos(V>) -C*fccos(6»)} 



k.b(p) [u(q).b(k)]{C*fccos(^) + C*fccos(6>)} ) 



(B.35) 



(B.36) 



We now take the limit k, q — > oo, keeping ip and p constant. In this limit <f> — > 0. Making use of the relationships 
9 = tt — ip — cj), k/ svn9 = p/sin</>, and taking the limit <fi — > 0, we get for Ya — ► k.b(p) [u(k).b(q)]{— Cpsin 2 (-0)} + 
k.b(p) [u(q).b(k)]{Cpsin 2 (V')}. However, in this limit Ya should tend to zero in order to satisfy locality. From this 



constraint we get C = 0. Putting the value of C in Eq. (B.34) gives — 



Therefore we get Ya = 0. It follows that the mode-to-mode transfer 7H bb = $ bb . 

In the above proof we used rotational invariance, galilean invariance, finiteness of mode-to-mode transfer, the 
symmetry of Ya with respect to k, p, q, u(k), u(p), u(q), b(k), b(p), b(q), and the assumption that the 
circulating transfer should vanish if two of the wavenumbers approach infinity. This proof is valid only upto linear 
orders, as explained earlier in this appendix. A completely general proof without taking recourse to the linearity 
assumption is lacking at present. 
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APPENDIX C: DERIVATION OF THE CIRCULATING TRANSFER BETWEEN VELOCITY MODES 

AND MAGNETIC MODES IN A TRIAD 



In Section VII C we had shown that the energy transfer from u(p) to b(k) in the triad (k, p, q) can be expressed 
as ^. 6tI (k|p|q) = 5# b "(k|p|q) + ^ & "(k|p|q) is given by Eq. ( f76| ) of Section |VII C was interpreted as the effective 
transfer from u(p) to b(k) and was interpreted as a circulating transfer between the velocity and magnetic 
modes belonging to a triad. We will follow the approach of Appendix A and Appendix B to determine Za in this 
Appendix. We will impose the following conditions on Z&. It should not violate the fundamentally triad nature of 
the interactions, it should be finite, rotationally invariant, galilean invariant, and should vanish if k/p — > oo. 

The circulating transfer Z/\ is a scalar with a dependence on atmost k, p, q, u(k), u(p), u(q), b(k), b(p), b(q). 



We know from Appendix B that such a function can be expressed in the form given by Eqs. (B.l)-( p.5| ) of Appendix 
B. The interactions between the velocity and the magnetic modes are fundamentally three-mode, i.e., each term in 
the expression for combined energy transfer to a velocity or a magnetic mode depends on all the three wavenumbers 
(k, p, q). the combined energy transfer to any of the velocity or the magnetic modes in the triad. Za should not 



violate this fundamental nature of the interactions. Hence we should drop all those terms from Eqs. (B.1)-(B.5) which 
are independent of atleast one of the modes and then we get Eqs. ( |B.6| )-( |B~^ ) for Z^. The combined energy transfer 
to velocity and magnetic modes in the triad are zero if {b(k), b(p) = 0}, {b(k), b(q) = 0}, or {b(p), b(q) = 
0}. The terms which are non- vanishing if one of the above set of Fourier coefficients are zero will be dropped from 



Eqs. (B L 6)-(B_ L 9), giving Eqs. (|B.10| )-( B.12 ) In Appendix B also we had imposed identical constraints on Eqs. ( B . 1 ) 



(B.5) to get Eqs. ( B.10 )-( B.12 ). The dimensionless coefficients in these equations are k, p, q, u(k), u(p), u(q), 
b(k), b(p), b(q) dependent and we will now attempt to determine these. 



In Appendix B we had considered the consequences of finiteness and galilean invariance on Eqs. ( B.10 )-( B.12 ) 
of Appendix B. We will briefly summarise the conclusions now. To ensure the finiteness and galilean invariance 



of Eqs. ( B.1C )-( B.12 ) the dimensionless coefficients in these equations should be independent of the magnitudes 
of the wavenumbers and the Fourier coefficients. Then each coefficient in the equations can be written as a — 



a (o) _|_ a (l) _|_ a (2) + a (3) + a (4) wne re a (0) is a constant, a (2) , a (3) defined in Eqs. jB.14j )-( pjq ) of Appendix 
B. are linear in terms of the kind k.p, k.u(p), u(k).u(p) respectively, and ay^ is nonlinear in these terms, a^- 2 ', 
and c^ 3 ) are not galilean invariant and would therefore lead to the non-invariance of Eqs. ( jB.14D -( [Bl6| ) of Appendix 
B. Hence, and should not be included in a. In Appendix B we had also dropped a^ 4 -* . Therefore, we get 
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a = + aS 1 '. Following Appendix B we will now write Z A [given by Eqs. (B.1C)-(B.12) in Appendix B] 



where 



Z% = Z& + Z& (C.l) 



Z 2 b A (k|p|q) = fe([k.u(q)][b(k).b(p)]{a<°) +4 1)3 ^ + + 

[k.u(p)] [b(k).b(q)] {/?(») + % + /#> ^ + /f 5 ^ } + 
[p.u(k)][b(p).b(q)]{ 7 (°) + 7 f )k f + 7^ +7f ) ^} ) (C2) 



i? e ([k.b(q)][u(k).b(p)]{^°)- 


H^ 5 


kp 


h4 1} 


kq j 
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pqi 

P9 J 


[k.b(q)][u(p).b(k)]{^ )- 
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kp 
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kq 




pqi 

P9 J 


[k.b(p)][u(k).b(q)]{ M (°)H 
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kp 




k^q 

kq 




pqi 

p<? ^ 


[k.b(p)][u(q).b(k)]{z/ ) - 




k^p 

kp 




k^q 

kq 




pqi 

P9 J 


[p.b(k)][u(p).b(q)]{^ 


-4P 


hj? _ 

kp 


^ 2 1} 


k^q 

kq 




pqi 

P9 



[p.b(k)] [u(q) .b(p)] {*«» + 5* + ^ + ^ } ) (C3) 



The quantities Z A " = + Zg A represent the energy circulating from u(p) — » b(k) — ► u(q) — > b(p)... — > u(p) [see 
Fig. 9 in Section VII C| |. We will denote the circulating transfer from u(p) — > b(k) by Z A u (k|p|q), from b(p) — * u(k) 



by Z A b (k|p|q) and so on for the circulating transfer between other modes. We take Z h £ = Z2& + bu Z^\ in Eqs. ( p.2| )- 
( |C. 3D to be the circulating transfer from u(p) — > b(k), i.e., Z A "(k|p|q) = Z A - Then by symmetry the circulating 
transfer from u(q) to b(k) can be written as 

Z£(k|q|p) = Z 2 b A (k|q|p) + Z 3 6 A (k|q|p) (C.4) 



Z 2 b A (k|q|p) = fe([k.u(p)][b(k).b(q)]{a(°) +ai ^+a 2 ^ + a 3 ^} + 
[k.u(q)][b(k).b(p)]{/?(°> +/3 1 k ^ +/3 2 ^ +/3 3 ^} + 
[ q .u(k)][b(p).b(q)]{ 7 ( )+ 7l ^+ 72 k ^ + 73 ^} ) (C.5) 

Z 3 fc A (k|q|p) = Re( [k.b(p)][u(k).b(q)]{<5(°) +S^ + 6 2 ±z + 5 3 ^f} + 
[k.b(p)][u(q).b(k)]{ ?7 ( ) + m ^ +m ^ + ^£3} + 
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[k.b(q)][u(q).b(k)]{ M (°) + n** + + M3^} + 

(k.b(q)][u(k).b(q)]{i/(°) + ^ + ^ + i***} + 
[q.b(k)][u(q).b(p)]{0(°) + ^ + fa** + h^f} + 
[q.b(k)][u(p).b(q)]{x(°)+Xi^ + X2^ + X3^} ) (C.6) 

We have dropped the superscript '1' in the above equations. Z^(k|p|q) and Z^ b (p|k|q) are physically the same 
transfers, but the former represents the transfer from u(p) to b(k) and the latter represents the same transfer from 
b(k) to u(p) The former should therefore be equal and opposite in sign to the latter, i.e., 

Z£(k|p|q) + Z£ b (p|k|q) = (C.7) 



(p|k|q) is thus given by the negative of Eqs. (C.1)-( |C.3| ). Using the expression for Z% (p|k|q) we can write 



Z£ & (k|p|q) to be 

Z£ b (k|p|q) = Z 2 lb (k|p|q) + Z 3 t(k|p|q) (C.8) 

Z 2 t(k|p|q) = Re( - [k.u( P )][b(k).b(q)]{ 7 (°) + 7l i|£ + 72 ^ + 73 ^} 
+ [k.u(q)][b(k).b(p)]{a(°) + a x % + a 2 ^f + a 3 ^} 
[p.u(k)][b(p).b(q)]{/3(°) + +/3 2 ^ + f 3 3 ±f}) 

(C.9) 
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kp 







(C.IO) 



and Z^ (k|q|p) is given by 
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ZX 6 (k|q|p) = Z 2 t(k|q|p) + Z3i(k|q|p) (C.ll) 

Z«i(k|q|p) = Re( [k.u(q)][b(k).b(p)]{ 7 W +71^+72^ +73^} 
+ [k.u(p)][b(k).b(q)]{a(°) + ai ±2+a 2 ^+ a 3 ^} 

(C.12) 



^3A(k|q|p) = 
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-Xl 


kq 
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kq 
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kq 
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z P9 


^3^}) 



(C.13) 



(C.14) 



The sum of the circulating transfer to b(k) from u(p) and from u(q) should be zero, i.e., 

Z£(k|p|q) + Z£(k|q|p)=0 

and similarly the sum of the circulating transfer to u(k) from b(p) and from b(q) should also be equal to zero, i.e., 

^ b (k|p|q) + Z£ b (k| q |p) = (C.15) 



Replacing Z^ l (k|p|q) and Z^ l (k|q|p) from Eqs. jcl|)-(|C^) and Eqs. (^J)-(^|) into Eq. ( |C.14| ) we get the 
following relationships between the constants: a, ft etc. + ft ) — 0, ot\ + ft = 0, a 2 + ft = 0, a 3 + ft — 

0, 7i - 72 = 0, <?(°) + = 0, <5i + fi 2 = 0, 6 2 + Mi = 0, 6 3 + M3 = 0, ?/ 0) + = 0,r}i + u 2 = 0, m +Pi = 0, 773 + 2/3 = 
0, (pW - x (0) = Q, 4>i - X2 = 0, 02 - Xi = 0, 03 - X3 = 0. After replacing ZX b (k|p|q) and ZX b (k|q|p) from Eqs. jCj|)- 



(|0§) and Eqs. ( p.ll| )-( |013| ) into Eq. dC.15| ) we get the following relationships: a 2 - a 3 = 0, ft ) +7 (0) = 0, ft + 71 
0, ft + 73 = 0, ft + 72 = 0, - =0,Si-rn = 0, <5 2 - % = 0, 5 3 - m = 0, M (0) + (o) = 0, f n + 0i = 0, /i 2 + 03 
0, Ai 3 + 02 = 0, t/(°) + X (0) = 0, 14, + xs = 0, M3 + 02 = 0. 
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Let us now consider the following geometry: the three wavenumbers in the triad form an isosceles triangle, and we 
choose Fourier coefficients at the two wavenumbers of equal lengths to be perpendicular to k, p, q. We had considered 
this kind of geometry in Appendix A and Appendix B. We had argued that because of the symmetry between the 
two wavenumbers of equal length, the circulating transfer should vanish for this geometry. This constraint had given 
us an additional set of relationships between the constants in Ya . Z h £ should also be equal to zero for this geometry. 



Since Z h £ and Ya have the same functional form, given by Eq. (B.22)-(B.23) of Appendix B and Eq. ( p.l[ )-(C.3) of 
this appendix respectively, hence imposition of this constraint on Z b £ would give the same set of relationships, i.e., 
7 (°) = o, 7l + 72 = 0, 73 = 0, = 0, 2 + 03 = 0, 0i = 0, x (0) = 0, xi + X3 - 0, X 2 = 0, a<°> = 0, ai = 0, a 2 + a 3 = 
0,£(°) = 0, 6 1 = 0, S 2 + 5 3 = 0, t?(°) = 0,77a = 0,m+V3 =0,pM = 0, & = 0, ft + /3 3 = 0, /x (0) = 0, M 2 = 0,1* +M3 = 
0, = 0, v 2 = 0, v x + v z = 0. 



The above relationships between the constants in Eqs. (G.2)-(C3) are identical to the relationships between the 



constants in Eqs. (B.22)-(B.23) of Appendix B. We know from Appendix B that these relationships can be expressed 



more clearly as follows: { a 



(0) 



-pW = 7 (o) = o}, {Si 



o) 



t (0) = ^(O) = ^(0) = _ i/ (o) 



AO) 



0}, {ai 



-ih 



7 3 = 0}, {a 2 = -p\ = 71 = 72 = a 3 = -/3 3 = 0}, {Si = -/j, 2 = = X3 = -v% = Vi = 0}, {<5 2 = -fii = 4>i = X2 = 
—^3 = 773 = —^3 = /x 3 = —02 = — Xi = = — 772}- There is just one undetermined constant which we will denote by 
'C, i.e., {82 = — Mi = •■■ = C}. Expressing Z b £ [Eqs. ( |C.2[ )-( |C^3| )] in terms of this single constant, we get 

Z 2A = (C.16) 

Z 3A = Re( [k.b(q)][u(k).b(p)]{C^ - C^} + [k.b(q)][u(p).b(k)]{-C^ + C ^} + 
[k.b(p)][u(k).b(q)]{-C^ + C^} + [k.b(p)][u(q).b(k)]{C^ + -C^} + 
[p.b(k)][u(p).b(q)]{C^ - C^} + [p.b(k)][u(q).b(p)]{-C^ + C^} ) 

(C.17) 



which is the same as the expression for for Ya in Eqs. (B.33)-(B.34) of Appendix B. 

If k/p and q/p — * 00, while p is kept constant the circulating transfer Za should tend to zero. This constraint was 
also applied on the circulating transfer between magnetic modes, Ya in Appendix B. We had shown that the constant 



C = for this constraint to be satisfied. Putting C = in Eqs. ( p.!6| )-( pT7| ) we get Z b £ = Z b \ + Z\\ equal to zero. 

The circulating transfer between the velocity and the magnetic modes in the triad is equal to zero and the mode- 
to-mode transfer ^. b "(k|p|q) = ^ h "(k|p|q). 
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The derivation of Z b £ , however took the coefficients to be linear. A completely general derivation for Z b £ will 
require imposition of additional symmetry and invariance or some other constraints, which we have ignored here. 



45 



APPENDIX D: EQUIVALENCE OF KRAICHNAN'S FORMULA AND OUR FORMULA FOR THE 

KINETIC ENERGY FLUX 



In this appendix we will show that the flux formula (30) based on the effective mode-to-mode transfer is equivalent 
to Eq. (^) derived by Kraichnan. 

Kraichnan showed that the kinetic energy flux in terms of the combined energy transfer S uu given by 

A 1 

K<(K) = - ]T £-S""(k|p,q) (D.l) 

k <K P 

Our formula for flux in terms of the effective mode-to-mode transfer fi uw is given by 

k>(k) = - E E ^"( k ip|q)- (d-2) 

|k|<if |p|>if 



The equation ( D.l ) can be written in terms of as follows: 

A 1 

K<(K) = E E o (^""(klPlq) + r"(k|q|p)) (D.3) 



|k|<K P 



which can be further split up as 



E ^(r U (k|p|q)+r"(k|q|p)) 
|k|<K|p|,|q|<if 

+ E E ^(r U (k|p|q)+r"(k|q|p)) 
|k|<K |p|,|q|>if 

+ E E ^(r"(k| P |q)+r"(k|q|p)) 

|k|<K |p|<K,|q|>iT 

+ E E 5(r"(k|p|q)+r"(k|q|p)) (D.4) 

|k|<K \p\>K,\q\<K 

These transfers are shown in Fig. D.l. Making use of symmetries between k, p,q in the summations we get the 
following relationships. Since p, q are interchangeable we get 

E E r"(k|p|q) = E E ^(k|q|p), (D.5) 
|k|<K |p|,|q|>K |k|<K|p|,|q|>if 

since k, p are interchangeable, 

E E r"(k|p|q) = E E ^"(Plklq), (D.6) 
|k|<K |p|<J<r,|q|>K |k|<K |p|<K,|q|>iC 

4G 



since k, q are interchangeable, 

E E ^"( k iq|p)= E E r u ^\k\ P ). (d.7) 

|k|<K \q\<K,\p\>K |k|<K |q|<A',|p|>A' 

The sum X]|k|<K X^|<if \^\>k ^-j^" 14 (k|p|q) involves transfers from k — + p and also from p — > k, which are equal and 
opposite. Therefore 



E E ^"(k| P |q) = 0. (D.8) 



|k|<K |p|<K,|q|>K 



Similarly, 



E E ^""(k|q|p) = 0, (D.9) 



2 

|k|<K |p|>K,|q|<K 



and, 



E E ^ U "(k|q| P )=0, (D.10) 

|k|<K|p|<K,|q|<iC 

In other words all terms except both legs of Fig. D.l, p — > k legs of Fig. D.l, and q — > k legs of Fig. D.l survive. 
We use these results to rewrite II"> (K) as 

A A 

2 7 



n:<(#)=E E r u (k| P |q)+E E ^ MU (k|q| P ) 

|k|<K|p|,|q|>K |k|<K|p|<if,|q|>K 

+ E E ^""(k|p|q) (D.ll) 



T 

|k|<K |q|<K,|p|>_R- 



Using symmetry between p and q we get, 



A A 



E E ^""(k|q|p) = E E ^ Wl (k| P |q) (D.12) 

Hence, we can write, 



|k|<K |p|<K,|q|>K |k|<K |q|<K,|p|>iC 



n ">=E E r"(k| P |q)+E E r"(k|p|q) 
|k|<K |p|,|q|>K |k|<K |q|<K,|p|<JC 

A 

= E E r"(k| P |q) (D.13) 
|k|<K \p\>K 

This is the same as of Eq. (p8[). Thus, we have shown that Kraichnan's formula for flux is the same as our 
formula based on the mode-to-mode transfer. 
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FIGURE CAPTIONS 

Fig. 1 two types of triads involved in transfer between shell m and shell n. Mode q could be either inside a shell or 
could be located outside the shells. 



Fig. 2 The mode-to-mode energy transfers T/t uu 's sought to be determined. 

Fig. 3 Mode-to-mode transfer can be expressed as a sum of circulating transfer and effective mode-to-mode transfer. 

Fig. 4 The circulating transfer and mode-to-mode transfer from shell n to shell m with the arrows pointing towards 
the direction of energy transfer. 

Fig. 5 The circulating transfer Xa does not contribute to the energy flux out of the sphere K. 

Fig. 6 The combined energy transfers in a triad. 5" u (k|p, q): u(p) and u(q) — > u(k); S bu (q\k,p): u(k) and 
u(p) - b(q); S 6b (k|p,q): b(p) and b(q) - b(k). 

Fig. 7 The mode-to-mode energy transfers sought to be determined. 

Fig. 8 Energy transfer between a pair of magnetic modes in a triad can be expressed as a sum of circulating transfer 
and effective mode-to-mode transfer. 

Fig. 9 Energy transfer between a velocity and a magnetic mode in a triad can be expressed as a sum of circulating 
transfer and effective mode-to-mode transfer. 

Fig. 10 The various energy fluxes defined in this section are shown. The arrows indicate the directions of energy 
transfers corresponding to a positive flux. These fluxes are computed in numerical simulations and the results are 
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given in a subsequent paper. 



Fig. B.l In the wavenumber triad k, q — > oo >> p the circulating transfer, Ya — to avoid non-local transfer from 
very large wavenumbers to very small wavenumbers in comparison. 



Fig. D.l he various triads involved in the terms in Eq. ( D.4 ). Triad of type A does not contribute to energy flux. 
Triad of type C and D are equivalent and contribute the same amount. 
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FIG. 1. 
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FIG. 2. 



•52 



U(q) 




u(k) uCp) 



Fig. 3 



FIG. 3. 
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Fig, A 



FIG. 4. 
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Fig . 5 



FIG. 5. 
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FIG. 6. 
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FIG. 7. 
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FIG. 8. 
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FIG. 10. 
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Fig. B. 1 
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Fig. D. 1 
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